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ABSTRACT 

In this first paper in a series we present 1298 low-redshift {z < 0.2) optical spectra 
of 582 Type la supernovae (SNe la) observed from 1989 through 2008 as part of the 
Berkeley SN la Program (BSNIP). 584 spectra of 199 SNe la have well-calibrated light 
curves with measured distance moduli, and many of the spectra have been corrected 
for host-galaxy contamination. Most of the data were obtained using the Kast double 
spectrograph mounted on the Shane 3 m telescope at Lick Observatory and have a 
typical wavelength range of 3300-10,400 A, roughly twice as wide as spectra from 
most previously published datasets. We present our observing and reduction proce- 
dures, and we describe the resulting SN Database (SNDB), which will be an online, 
public, searchable database containing all of our fully reduced spectra and compan- 
ion photometry. In addition, we discuss our spectral classification scheme (using the 
SuperNova IDentification code, SNID; Blondin & Tonry 2007), utilizing our newly 
constructed set of SNID spectral templates. These templates allow us to accurately 
classify our entire dataset, and by doing so we are able to reclassify a handful of ob- 
jects as bona fide SNe la and a few other objects as members of some of the peculiar 
SN la subtypes. In fact, our dataset includes spectra of nearly 90 spectroscopically 
peculiar SNe la. We also present spectroscopic host-galaxy redshifts of some SNe la 
where these values were previously unknown. The sheer size of the BSNIP dataset and 
the consistency of our observation and reduction methods makes this sample unique 
among all other published SN la datasets and is complementary in many ways to the 
large, low-redshift SN la spectra presented by Matheson et al. 2008 and Blondin et al. 
(in preparation). In other BSNIP papers in this series, we will use these data to ex- 
amine the relationships between spectroscopic characteristics and various observables 
such as photometric and host-galaxy properties. 

Key words: cosmology: observations - distance scale - supernovae: general - surveys 



1 INTRODUCTION 

Supernovae (SNe) have been integral to our understanding 
of the cosmos throughout the history of astronomy - from 
demonstrating that the sky was not unchanging beyond the 



lunar sphere (Brahe 15731 to the discovery of the acceler- 



ation of the expansion of the Universe (Riess et al. 1998 



Perlmutter et al.||1999"l ). Type la supernovae (SNe la) have 



been particularly useful in recent years as a way to accu- 



rately measure cosmological parameters 


Astier et al. 


2006 


Riess et al.|2007 Wood-Vasey et al. 


2007 


Hicken et al. 


2009 


Kessler et al.|2009||Amanullah et al. 


2010 


Conley et al. 


2011 


Sullivan et al.||2011 


Suzuki et al.||2011 


I. Broadly speaking 



SNe la are the result of thermonuclear explosions of car- 



bon/oxygen white dwarfs (WDs) (e.g.. 


Hoyle & Fowler] 1960| 


Colgate & McKee|1969||Nomoto et al. 


1984| see|Hillebrandt| 


& Niemeyer|2000 for a review). However, we still lack a de- 



tailed understanding of the progenitor systems and explo- 
sion mechanisms, as well as how differences in initial condi- 
tions create the variance in observed properties of SNe la. To 
solve these problems, and others, detailed and self-consistent 
observations of many hundreds of SNe la are required. 

The cosmological application of SNe la as precise dis- 
tance indicators relies on being able to standardise their lu- 
minosity. |Phillips| ( |1993[ ) showed that light-curve decline is 
well correlated with luminosity at peak brightness for most 
SNe la, the so-called "Phillips relation." Additionally, people 
have searched for a "second parameter" in SN observations 
which would make our measurements of the distances to 



SNe la even more precise. Bailey et al. (20091 and Blondin 
|et al.| ( |2011[ ) have decreased the scatter in residuals to the 
Hubble diagram with the help of spectroscopic data. |Wang| 



et al. (20091 obtained an additional improvement by sep- 
arating their sample of SNe la into two groups based on 
the ejecta velocity near maximum brightness; they suggested 
different reddening laws for these two samples. Building on 
this work, Foley & Kasen (20111 found that the intrinsic 



maximum-light colour of SNe la depends on their ejecta ve- 
locity. After accounting for this colour difference, the scatter 
in Hubble-diagram residuals is decreased from 0.19 mag to 
0.13 mag for a subset of SNe la. This particular conclusion 
was possible only with a large set of spectroscopically ob- 
served objects, with many of the spectra coming from the 
sample described in this Paper (see also Wang et al.|2009 l. 

Until now there have been several statistical samples 
of low-redshift SN la photometry (e.g., Hamuy et al.||1996 



Riess et al.||T999| [Jhaet al.|2006| [Hicken et al.l2009| IGarie 
shalinga m et al.|20To ContrerasetaLpOlOI IStritzinger et al. 
2011|), but only one large sample of low-redshift SN la spec- 
tra ( Matheson et al.|2008 1 . Until the publication of over 432 
spectra of 32 SNe la by Matheson et al. ( 2008'), large sam- 



ples of SN la spectra were typically constructed by combin- 
ing datasets published for individual objects, usually from 
many different groups. 

The Berkeley Supernova la Program (BSNIP) is a 
large-scale effort to measure the properties of low-redshift 
{z < 0.2) SNe la, focusing on optical spectroscopy and pho- 



tometry (see Ganeshalingam et al. 2010 for the compan- 



ion photometry paper to much of the spectroscopic sam- 
ple presented here). One aspect of our strategy for the last 
two decades has been to observe as many SNe la as possi- 
ble in order to dramatically increase the number of objects 
with spectroscopic data. We have also attempted to obtain 



BSNIP I: SN la Spectra 3 



good temporal spectral coverage of peculiar objects as well 
as objects which were being observed photometrically by our 
group. In addition, we strove to spectroscopically classify all 
SNe discovered by the 0.76 m Katzman Automatic Imaging 
Telescope (KAIT; Filippenko et al.|2001 1. By observing and 
reducing our spectra in a consistent manner, we avoid many 
of the systematic differences found in previous samples con- 
structed from data obtained by various groups. 

In this paper we present the low-redshift SN la spectral 
dataset. This sample consists of a total of 1298 spectra of 
582 SNe la observed from 1989 through the end of 2008. A 
subset of the SNe, along with information about their host 
galaxies, is presented in Table [l] (the full set is available 
online - see Supporting Information) . Information regarding 
some of the SN la spectra in the dataset is listed in Table [2] 
(and again the full set is available online - see Supporting 
Information). Many spectra presented in this paper have 



complementary light curves from Hamuy et al. ( 19961, Riess 



et al. (19991, Jha et aL (20061, Hicken et al. (20091, and 



Ganeshalingam et al. (20101, which have all been compiled 
and fit by Ganeshalingam et al. (in prep.). Other spectra 
have complementary unfiltered light curves given by Wang 
et al. (in prep.). 

In this paper, we describe our observations and data- 
reduction procedure in Sections [2] and [3] respectively. We 
present our methods of data management and storage in 
Section [4] and our spectral classification scheme in Section [5] 
The sample of objects and spectra is described in Section^ 
and there we also show our fully reduced spectra as well eis 
(for the objects with multi-band SN and galaxy photome- 
try) galaxy-subtracted spectra. Reclassifications of a handful 
of SNe, as well as previously unknown spectroscopic host- 
galaxy redshifts, are also given. We discuss our conclusions 
in Section]?] Future BSNIP papers will examine the correla- 
tions between spectroscopic properties and other observables 
(such as photometry and host-galaxy properties). 
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Figure 1. The percentages of spectra presented in this work ob- 
tained using various instruments. The UV Schmidt spectrograph 
( [Miller fc Stone|1987| l was on the Shane 3 m telescope at Lick Ob- 
servatory until early 1992 when it was replaced by the Kast dou- 
ble spectr ograph ([Miller fc Stone|19 93). LRI S l[Oke et al.|1995] , 
DEIMOS l [Faber et al.|2003^ , andESI (Shcinis et al.|2002| were 
all used with one of the 10 m Keck telescopes. 



for each object (see Section 6.1 for further comparisons of 
the two spectral datasets). 

All observations of our scheduled time were performed 
by members of the BSNIP group. PI Filippenko was present 
for an astounding 254 nights. Occasionally, as the result of 
a swap of time or for a particularly interesting object, an 
observer exterior to the BSNIP group would observe for our 
team. This sometimes resulted in slight variations in instru- 
ment configurations (such as a smaller wavelength range, for 
example). As mentioned above, the bulk of our data were ob- 
tained at the Lick and Keck Observatories where our average 
seeing was slightly greater than 2" and slightly greater than 
1", respectively. 

2.1 Individual Instruments 

2.1.1 UV Schmidt on the Shane 3 m 

The UV Schmidt spectrograph contained a Texas Instru- 
ments 800 X 300 pixel charge-coupled device (CCD) and our 
setup used a slit that was 2-3" wide. The average resolution 
of our spectra from this instrument was ~12 A. 



2 OBSERVATIONS 

Over the past two decades, our group has had access to sev- 
eral different telescopes and spectrographs for the purpose 
of observing SNe. The main facility for this study was the 
Shane 3 m telescope at Lick Observatory. During this pe- 
riod, the Shane telescope has had two low-resolution spec- 
trographs: the UV Schmidt spectrograph until early 1992 



(Miller & Stone 19871 and the Kast double spectrograph 
since then ( Miller fc Stone|1993 |. We also obtained a hand- 
ful of spectra using the Stover spectrograph mounted on the 
Nickel 1 m telescope also at Lick. 

We have supplemented our Lick Observatory sample 
with spectra obtained at the Keck Observatory. When con- 
ditions were not acceptable for our faint, primary targets 
(typically in twilight, or during times of bad seeing or 
cloudy weather), we would use one of the 10 m Keck tele- 
scopes to obtain spectra of our relatively bright (typically 
7? < 18 mag), nearby SN targets. We also obtained many 
late-time spectra with the Keck telescopes. We have used 



the Low Resolution Imaging Spectrometer (LRIS; Oke et al 



1995 



both before and after the addition of the blue arm), 
the DEep Imagin g Multi-Object Spectrograph (DEIMOS; 
Faber et al. 20031, and the Echelle Spectrograph and Im- 
ager (ESI; Sheinis et al.|[2002 |. The percentages of spectra 
presented in this work obtained using these instruments are 
displayed in Figure [T] 

All of these telescopes were classically scheduled. We 
would typically have 1 night every two weeks on the Shane 
telescope (near first and last quarter moon) throughout the 
year and 4-10 nights per year with the Keck telescopes (typ- 
ically 1-2 nights near new moon in a given lunation). Re- 
cently, we have been allotted a third night per lunar cycle 
on the Shane telescope near new moon. Taking into account 
weather and instrument problems, our coverage of any given 
object is typically about one spectrum every two weeks. The 
telescope scheduling and observing method are very different 
from those of'Matheso n et aL] ( |2008| l and Blondin et al. (in 
prep.), who observed fewer SNe la but with a higher cadence 



2.1.2 Kast on the Shane 3 m 

Until September 2008, the Kast double spectrograph used 
two Reticon 1200 x 400 pixel CCDs with 27 pixels and a 
spatial scale of 0.78" pixeP^, with one CCD in each of the 
red and blue arms of the spectrograph. Currently, the blue 
arm of Kast uses a Fairchild 2048 x 2048 pixel device with 
15 fim pixels, which corresponds to 0.43" pixel"^. For our 
typical setup, we would observe with a 300/7500 grating for 
the red side, a 600/4310 grism for the blue side, and a D55 
dichroic. This results in a wavelength range of 3300-10,400 A 
with overlap between the two arms of 5200-5500 A. With 
our typical slit of 2", we achieve a resolution of ~11 A and 
~6 A on the red and blue sides, respectively. 



2.1.3 Stover on the Nickel 1 m 

The Stover spectrograph contains a Reticon 400 x 1200 pixel 
CCD and 27 fim pixels with a spatial scale of 2" pixel^^. 
Our setup used a 2.9" wide slit with the 600/4820 grism. 
This yielded an average resolution of ~7 A. 



2.1.4 LRIS on the Keck 10 m 

When most of our dataset was obtained, the LRIS spec- 
trograph used a Tektronix 2048 x 2048 pixel CCD with 
21 fim pixels and a spatial scale of 0.211" pixel"'^ for the 
red arm and two 2048 x 4096 pixel Marconi E2V CCDs with 
15 fim pixels and a spatial scale of 0.135" pixel"^ for the 
blue arm. LRIS operated with only the red arm until 2000. 
The original blue-side CCD, used from 2000 to 2002, was an 
engineering-grade SITe 2048 x 2048 pixel CCD. Our typical 
setup would use the 400/8500 grating for the red side, either 
the 400/3400 or 600/4000 grism for the blue side, and the 
D56 dichroic, resulting in a wavelength range of 3050-9200 A 
and 3200-9200 A for the respective grisms. There was typi- 
cally an overlap region of 5400-5800 A and 5400-5700 A for 
the 400/3400 and 600/4000 grisms, respectively. With our 
typical 1" slit, this setup yields resolutions of ~7 A for the 
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red side, and either ~6.5 A or ~4.5 A for the 400/3400 and 
600/4000 grisms, respectively, for the blue side. 

2.1.5 DEIMOS on the Keck 10 m 

The DEIMOS spectrograph uses a 2 x 4 mosaic of 2048 x 4096 
pixel CCDs with 15 pixels and a spatial scale of 0.1185" 
pixel"^ for a total detector array of 8192 x 8192 pixels. 
Our typical setup would use the 600/7500 grating with a 
GG455 order-blocking filter, resulting in a wavelength range 
of 4500-9000 A. We would generally use a 1.1" slit which, 
along with our typical setup, would result in a resolution 
of ~3 A. Occasionally we would use the 1200/7500 grating 
instead, yielding a wavelength range of 4800-7400 A and a 
resolution of ~1.5 A. The slit was tilted slightly to provide 



better sky subtraction (see Section 3.1.4 for details) 



2.1.6 ESI on the Keck 10 m 

The ESI spectrograph has an MIT-LL 2048x4096 pixel CCD 
with 15 /im pixels and an average spatial scale of 0.1542" 
pixel" ^, with the redder orders having a larger spatial scale 
than the bluer orders. Our observations were typically per- 
formed in the echellette mode with a 1" slit in 2 x 1 bin- 
ning mode (spatial x spectral). This resulted in a resolution 
of 22 km s~^ across the entire wavelength range of 3900- 
11000 A. 




RC 
19.0% 



TM 
10.5% 



TS 
12.1% 



Figure 2. The percentages of spectra presented in this work 
reduced by various coauthors. JS, Jeffrey M. Silverman; RF, 
Ryan J. Foley; RC, Ryan Chornock; TS, Thea N. Steele; TM, 
Thomas Matheson; DL, Douglas C. Leonard; AB, Aaron J. Barth; 
O, Other. 



of our methods can be found in 



et ah (20011, Foley et al. (20031, Foley et al. (20071 



but the general method is the same. Previous descriptions 

Matheson et af] ([2000l ), p| 
and 

Matheson et al. (|2008|), but the discussion below supersedes 



previous descriptions. 



2.2 Standard Observing Procedure 



3.1 Calibration 



Unless there was a hardware malfunction, we would observe 
several dome flats at the beginning of each night (and oc- 
casionally at the end). We would also observe emission-line 
calibration lamps ( "arcs" ) at both the beginning of the night 
and often at the position of each object. Our final calibra- 
tions relied on observing standard stars throughout the night 
at a variety of airmasses. The goal was to obtain at least 
one standard star (in the case of single-beam spectrographs; 
both blue and red standard stars for double-beam spectro- 
graphs) at an airmass near 1.0 and at least one at an air- 
mass comparable to or higher than the highest airmass of 
any SN observed during that night. The standard stars were 



typically from the catalogs of Oke & Gunn ( 1983 1 and Oke 



(19901, with the cool, metal-poor subdwarfs and hot subd- 



warfs calibrating the red and blue sides, respectively. 

Most observations were made at the parallactic angle to 



reduce differential light loss (Filippenko 1982) . Exceptions 
were usually at an airmass < 1.2 or when the slit was po- 
sitioned at a specific angle to include a second object (the 
host-galaxy nucleus, a trace star, a second SN, etc.). In Au- 
gust 2007, LRIS was retrofitted with an atmospheric disper- 
sion corrector (ADC; |Phillips et al.|[2006| . With the ADC, 
differential light loss is substantially reduced regardless of 
position angle, even at high airmass. 



3 DATA REDUCTION 

All data were reduced in a similar, consistent manner by only 
a handful of people. The percentages of spectra presented 
in this work reduced by various coauthors are displayed in 
b igure [2] There are slight differences for each instrument. 



Despite differences between instruments, the general pro- 
cedure for transforming raw, two-dimensional spectrograms 
into fully reduced, wavelength and flux calibrated, one- 
dimensional spectra is similar for all of our data. We will 
discuss differences in the procedure for the various instru- 
ments below. The general prescription is as follows. 

(i) Correct for bias using an overscan region and trim the 
two-dimensional images to contain only the region with sky 
data. Our data do not typically show a bias pattern and do 
not have large dark currents. Therefore, we do not subtract 
bias frames, which would increase noise. 

(ii) Combine and normalise ffat-field exposures. We pay 
particular attention to masking emission lines from the fiat- 
field lamps and absorption features from the air between 
the fiat-field screen/dome and the detector. The normalizing 
function is generally a low-order spline. 

(iii) Correct pixel-to-pixel variations in our spectra using 
our fiat-field exposures. 

(iv) Extract the one-dimensional spectra. We use local 
background subtraction, attempting to remove as much 
host-galaxy contamination as possible. The spectra are typ- 
ically optimally extracted using the prescription of |Horne| 
(19861). 



However, prior to mid-1997, we did not use the optimal ex- 
traction for our Kast data but, rather, employed "standard" 
extractions (i.e., not optimally weighted). While this typi- 
cally had minimal impact on the signal-to- noise ratio (S/N) 
achieved (most SNe observed with Kast are quite bright, 
thus standard and optimal extractions yield similar noise 
levels in the final spectrum), one possible effect this has on 
our spectra from this era is that they may be spectrophoto- 
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metrically inaccurate at the ~5 per cent level. This is due 
to time- variable spatial focus variations that existed across 
the CCDs in the Kast spectrograph. By using the optimal 
extraction since then, we have significantly mitigated the 
effects of these variations in our data. 

(v) Calibrate the wavelength scale. Using arc-lamp spec- 
tra, we fit the wavelength scale with a low-order polynomial 
and linearise the wavelength solution. We then make small 
shifts in the wavelength scale to match the night-sky lines 
of each individual spectrum to a master sky template. 

(vi) Flux calibrate the spectra. We fit splines to the con- 
tinua of our standard-star spectra, producing a sensitivity 
function that maps CCD counts to flux at each wavelength. 
These sensitivity functions are then applied to each individ- 
ual SN spectrum. 

(vii) Correct for telluric absorption. Using our standard- 
star spectra, we interpolate over atmospheric absorption re- 
gions, providing an estimate of the atmospheric absorption 
at a particular time and airmass. Then, accounting for the 
differences in airmass, we apply these corrections to our 
spectra, allowing for slight wavelength shifts between the 
"A" and "B" telluric absorption bands. 

(viii) Remove cosmic rays (CRs) and make other minor 
cosmetic changes. In the remaining one-dimensional spectra 
there may be unphysical features due to CRs, chip gaps, or 
bad, uncorrected pixels. We interpolate over these features. 

(ix) Combine overlapping spectra. For instruments with 
both a red and blue side (or multiple orders in the case 
of ESI) we combine the spectra, scaling one side to match 
the other in the wavelength region where the spectra over- 
lap. For multiple, successive observations of the same object 
we combine the spectra to achieve the highest S/N in the 
resulting spectrum, weighting each spectrum appropriately 
(usually by exposure time). 

Through mid- 1997, we used our own VISTA and For- 
tran routines to complete all of the above steps. For about 
a year after that we used a combination of generic-purpose 
irafQ routines and our own Fortran routines for our spec- 
tral reductions. Since about mid-1998, we have performed 
our reductions using both generic-purpose IRAF routines 
and our own IDL scripts. Step 1 is achieved with either IRAF 
or IDL depending on the instrument. Steps 2-5 are gener- 
ally performed with IRAF, while steps 6-9 are performed in 
IDL. 

We consider the resulting spectra "fully reduced." How- 
ever, for a subsample of our spectra where we have multi- 
filtered host-galaxy photometry at the position of the SN 
and SN photometry near the time the spectrum was ob- 
tained, we can make additional corrections to obtain an ac- 
curate absolute fiux scale as well as account for host-galaxy 
contamination (see Sections |3.2| and |3.3|l . 



3.1.1 Kast on the Lick 3 m 

The Kast spectrograph has large amplitude, variable fring- 
ing on the red-side CCD. We observe red-side dome flats 
at the position of each object and we apply these fiats to 
each object individually. As the dome moves into place to 
take flats, we also obtain a red-side arc exposure. Using this 
arc spectrum we shift the wavelength solution derived from 
our afternoon arc exposures (which typically have more lines 
and are observed with a 0.5" slit, yielding higher-resolution 
lines) and apply those wavelength solutions to the appropri- 
ate SN observations. However, we still apply a small wave- 
length shift based on the night-sky lines later in the reduc- 
tion process. 



3.1.2 Stover on the Nickel 1 m 

The Stover spectrograph does not have the ability to rotate 
the slit with respect to the sky; thus, all spectra obtained 
with this instrument were observed with a fixed sky posi- 
tion angle of 0° . When observations were at relatively large 
airmasses (as they were for some of the spectra presented 
here) this caused their continuum shape to be unreliable. 
The spectra in our dataset from the Stover spectrograph 
have been previously published (Wells et al.||1994| |Li et al.] 
20011, and while strange spectrophotometric calibration is- 



sues when using this instrument with our setup and reduc- 
tion routines have been noted by [Leonard et ah] ( |2002[ ), [lI| 
et al. ( 2001 1 find no such problems. 



3.1.3 LRIS on the Keck 10 m 

The blue side of LRIS has two CCDs offset in the spatial 
direction (allowing a full spectrum to be on a single CCD). 
We typically position our objects on the slit so they will be 
centred on one CCD, ignoring the other CCD completely. 
However, in some observations circumstances dictated that 
objects be on the other CCD. Each CCD must be calibrated 
separately (different flat-field response functions, sensitiv- 
ity functions, etc.). The LRIS fiat-field lamp is not particu- 
larly hot, providing few photons at the bluest wavelengths 
of LRIS. We therefore mask this region in the fiat-field re- 
sponse, leaving the bluest portions uncorrected for pixel-to- 
pixel variations. 

When our dataset was obtained, the red-side CCD of 
the spectrograph had large fringes. We account for these 
fringes by applying dome fiats obtained during the afternoon 
or morning. We occasionally obtained internal flats at the 
position of an object, but we have found these to typically be 
worse for removing pixel-to-pixel variations than the nightly 
dome flats. However, there are rare instances where they 
were used instead of dome fiats. 



^ IRAF: The Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy (AURA), Inc., under cooperative agreement with the 
National Science Foundation (NSF). 



3.1.4 DEIMOS on the Keck 10 m 

The long slit for DEIMOS is slightly tilted, producing 
slightly different wavelengths for a pixel in a given column. 
This tilts the sky lines, providing additional sampling of 
the lines. Since our typical procedure is not adaptable to 
tilted sky lines (our background subtraction would produce 
dipoles for every sky line), we implement a modified version 
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of the DEEP2 DEIMOS pipelinfj^to rectify and background- 
subtract our spectra. The pipeline bias-corrects, flattens, 
traces the slit, and fits a two-dimensional wavelength so- 
lution to the slit by modeling the sky lines. This final step 
provides a wavelength for each pixel. The slit is then sky 
subtracted (in both dimensions) and rectified, producing a 
rectangular two-dimensional spectrum where each pixel in a 
given column has the same wavelength. From this point we 
resume our normal procedure, starting with extracting the 
spectrum (step 4). Since the spectrum has already been sky 
subtracted, we do not attempt any additional sky subtrac- 
tion (which would only increase the noise) unless the SN is 
severely contaminated by its host galaxy. 



in that column. We have found, however, that local sky sub- 
traction generally produces better results. 

On rare occasions, we have multiple dithered images of 
a single object. With these images we can (after proper scal- 
ing) subtract one from another to remove residual fringing 
and sky lines. We can also shift the spectra spatially and 
combine the two-dimensional spectra to increase the S/N of 
the object. This can produce better traces. 

For objects without a defined trace across the entire 
chip, we would create a trace function for the object either 
using the trace of a nearby object such as the host-galaxy 
nucleus or of a bright star (often an offset star) taken in the 
previous exposure at the same position as the SN. 



3.1.5 ESI on the Keck 10 m 



3.2 Spectrophotometry 



The CCD of ESI has several large defects which we mask 
before starting our reductions. These produce ~50 A gaps 
in our spectra, usually near 4500-4600 A. They can also 
affect our measurement of the trace, but this was rarely a 
problem for the low-redshift, relatively bright SNe presented 
in this paper. ESI observes 10 orders; we reduce each order 
individually and stitch the orders together at the end (step 
9), using our standard-star observations to determine the 
scaling and overlap regions. We weight the spectra in the 
overlap by their variance in each pixel before combining. 
For ESI we do not linearise the wavelength solution, but 
instead rebin to a common velocity interval, thus producing 
pixels of different sizes in wavelength space. 

3.1.6 Other Instruments 

In addition to the aforementioned instruments, our dataset 
contains a few spectra which were obtained by observers 
exterior to the BSNIP group at observatories aside from 
Lick and Keck. These data come from the Low Disper- 



sion Survey Spectrograph 3 (LDSS-3; Mulchaey & Gladders 
|2005[ ) mounted on the 6.5 m Clay Magellan II telescope, 
the R.C. spectrograph mounted on the Kitt Peak 4 m tele- 
scope, and the Double Spectrograph mounted on the Hale 
5 m telescope at Palomar Observatory | |Oke &: Gunn|1982| ). 

3.1.7 Additional Reduction Strategies 

Occasionally our standard procedures produce nonoptimal 
spectra. In these cases we augment our procedures to pro- 
duce higher-quality spectra. 

For some spectra we perform a CR cleaning of the two- 
dimensional spectra before extraction (step 4). This proce- 
dure is done in IRAF and detects pixels that have signifi- 
cantly more counts than their surrounding pixels, replacing 
them with the local median. Since this procedure has the 
potential to remove real spectral features, it is not automat- 
ically performed on every spectrum. 

We can obtain better sky subtraction on some spec- 
tra by performing a two-dimensional sky subtraction. This 
procedure fits each pixel in the spatial direction with a poly- 
nomial or spline function (usually constrained to the region 
near the SN position) and subtracts that fit from each pixel 



http : //astro .berkeley. edu/-cooper/deep/ spec2d/ 



Using our standard reduction procedure outlined in Sec- 
tion |3.1| the spectrophotometry of our data is usually quite 
accurate. However, there are many ways in which the spec- 
trophotometry of a spectrum may be corrupted. First, there 
are achromatic effects such as clouds that affect the abso- 
lute spectrophotometry. Absolute spectrophotometry is not 
necessary for many spectroscopic studies (although we will 
discuss absolute spectrophotometry in more detail in Sec- 
tion |3.2.2[ | . Accurate relative spectrophotometry, however, 
is important for such studies. There are many reasons why 
the relative photometry of a spectrum may be incorrect, but 
variable atmospheric absorption, nonparallactic slit angles 
leading to differential light losses ( Filippenko|1982 1 , and in- 
correct standard-star spectrophotometry can all contribute 
significant errors. As shown below, after rigorous testing we 
find that the relative spectrophotometry of the BSNIP data 
is accurate to ~0. 05-0.1 mag across most of the wavelengths 
covered by the spectra. 

Occasionally, events external to the normal operations 
of the telescope and data reduction can result in question- 
able spectrophotometry. Instrument failures (e.g., a broken 
shutter) or environmental effects (e.g., nearby wild fires) are 
the most troublesome. There is no clear way to fully correct 
the spectrophotometry in these cases. Using our detailed 
records as well as those of Lick Observatory, we have iden- 
tified several spectra where the spectrophotometry may be 
affected by these external factors and exclude them from any 
estimates of the fidelity of our spectrophotometry. Including 
spectra obtained with the Stover spectrograph, which does 
not have a rotator and so nearly all spectra were not ob- 
served at the parallactic angle, we have fiagged 88 spectra 
as having possibly troublesome spectrophotometry. 

3.2.1 Relative Spectrophotometry 

Two of the key attributes of the BSNIP sample are the 
large wavelength range and the consistent and thorough re- 
duction procedures. The spectra in the sample likely have 
similar systematic (and hopefully small) uncertainties. The 
large wavelength range makes the spectra ideal for compar- 
ing near-ultraviolet and near-infrared features in a single 
spectrum, but such investigations will be limited by the ac- 
curacy of our spectrophotometry. Since most of the spec- 
tra in our sample have corresponding BVRI light curves 
( [Ganeshalingam et al.||2010[ ), we can test the spectropho- 
tometry of a spectrum by comparing synthetic colours from 
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the spectrum to those of the hght curves at the time that 
the spectrum was obtained. In fact, this has previously been 
performed on some of the data presented herein (at a some- 
what less rigorous level) by Poznanski et al. (20021. 

For this test, we examine only the spectra of objects 
that have corresponding filtered light curves. To assure that 
our estimates of the SN colours from the photometry are 
accurate, we further limit the sample to spectra that have 
a light-curve point within 5 d of when the spectrum was 
taken. 

We use the light-curve fitter "Multi-colour Light Curve 
Shape" (MLCS2k2, |Jha et al.|[2007l ) to model the filtered 
light curves, allowing us to interpolate between data points. 
We fit each filter individually to provide the largest degree of 
flexibility in each, and all of the fits are inspected to ensure 
that a good fit is obtained. In cases where the MLCS2k2 
fit does not adequately reflect the data and the data are 
well sampled, we use a cubic spline with a Savitzky-Golay 
smoothing filter ( Savitzky fc Golay|1964 1 . 

We estimate the uncertainties in the model light curve 
by running a series of Monte Carlo simulations. For each 
data point, we randomly draw from a Gaussian distribution 
with mean given by the reported magnitude and a by the 
photometric uncertainty to produce a simulated data point. 
Each simulated light curve is refit. This process is repeated 
50 times and the scatter in the derived light curves is taken 
as the uncertainty in the model. This process is applied to 
objects with MLCS2k2 and spline fits. 

To determine the synthetic photometry from the spec- 
tra, we convolve each spectrum with the Bessell filter 



functions (Bessell 19901. We calibrate our photometry by 



measuring the spectrophotometry of the standard star 



BD-|-17°4708 (Oke & Gunn 19831 and applying zero-point 



offsets to match the standard photometry. We then apply 
these offsets to the synthetic photometry derived from the 
SN spectra. The Bessell filter functions have approximate 
wavelength ranges of 3700-5500, 4800-6900, 5600-8500, and 
7100-9100 A for B, V, R, and I, respectively. Most of our 
spectra fully cover the BVRI bands. 

There are several effects which may reduce the accu- 
racy of our spectrophotometry. By far, the most important 
is galaxy contamination. Although our reduction process re- 
moves as much galaxy hght as possible from a SN spectrum 



(see Section 3.1 1, some of our SN spectra are still contami- 
nated by galaxy light. The measured synthetic colours from 
galaxy-contaminated spectra will likely be vastly different 
from the SN colours even if our spectrophotometry is excel- 
lent. For spectra with multi-colour template images of the 
host galaxy and multi-colour light curves concurrent with 
the spectrum, we can correct for galaxy contamination to a 
large degree (see Section [3.3[ ). However, this correction relies 
on excellent relative spectrophotometry. 

We have selected a subsample of SNe that are relatively 
isolated from their host galaxy, so their spectra should have 
minimal galaxy contamination. These objects all have tem- 
plate images (taken after the SN had faded) that indicate 
minimal galaxy light. A sample of spectra of objects from 
this low-contamination sample of SNe is constructed to test 
the fidelity of our relative spectrophotometry. For this sam- 
ple, we require that the spectra have f < 30 d and that 
the spectrum was obtained at the parallactic angle or at an 
airmass ^ 1.2. We present the synthetic and photometric 



colours for the low-contamination sample in Figure [3] Al- 
though the number of spectra in this sample is limited, they 
span a large range of colour. 

We present a comparison of synthetic colours derived 
from our low-contamination and possibly contaminated 
spectra to those measured from light curves at the same 
epoch in Figure |4] An estimate of the uncertainty in the 
spectrophotometry can be made by examining the per de- 
gree of freedom (dof ) of the residual of the synthetic to pho- 
tometric colours. The uncertainty in the photometric colours 
is measured by examining the residuals of the photometry 
measurements near the epoch of the spectrum relative to 
the model. The uncertainty in the relative spectrophotom- 
etry is the uncertainty added to each point which causes 
the residuals of the synthetic to photometric colours to have 
X^/dof = 1. If x^/dof ^ 1 with only photometric uncertain- 
ties, then the spectrophotometry does not have uncertainties 
larger than the photometry itself. We present estimates of 
the uncertainties in Tabled 

For the low-contamination sample, the spectrophotome- 
try has a typical additional uncertainty of ^ 0.07 mag across 
the entire spectrum, with no additional uncertainty required 
for V — R and very little additional uncertainty (0.008 mag) 
required for R — I across a large range of colours. Our entire 
sample is only slightly worse, with the additional uncertainty 
inV — R being 0.055 mag. The accuracy of the flux calibra- 
tion for the low-contamination sample is difficult to assess 
since the uncertainties from the photometry are enough to 
account for the majority of the scatter (and the entire scatter 
for the wavelength region spanning from V to R). Nonethe- 
less, we can place limits on the accuracy based on the ad- 
ditional uncertainty required and the standard deviation. 
From this, we flnd that the low-contamination sample is ac- 
curate to 5.2-6.9 per cent, 0.0-5.8 per cent, 0.7-4.5 per cent, 
and 6.0-9.0 per cent for the wavelength regions spanning B 
to V, V to R, R to I, and B to I, respectively. For the sample 
of objects corrected for galaxy contamination, the additional 
errors are similar to those of the low-contamination sample 
(5.3-6.5 per cent, 3.9-4.8 per cent, 4.9-5.1 per cent, and 
4.5-9.2 per cent for the wavelength regions listed above), 
but lower than those for the entire sample (8.8 per cent, 
5.1-6.2 per cent, 7.3-8.9 per cent, and 12.7-15.6 per cent), 
indicating that the galaxy-contamination correction works 
well at least for broad-band colours. 

We have also split our sample by various spectral at- 
tributes. The spectrophotometry does not depend signifl- 
cantly on airmass. It does depend signiflcantly on S/N, but 
the spectrophotometry does not improve as S/N increases 
beyond S/N ~ 20 pixel" ^. The additional uncertainties also 
depend slightly on the individual who reduced the spectra. 
However, this trend may be the result of observation and 
reduction techniques slowly improving over time. 

We have also calculated the mean and standard devia- 
tions of the difference between the synthetic colours derived 
from our spectra and those measured from light curves for 
the various subsamples. All subsamples have a mean that 
is < 0.6 standard deviations from zero, with nearly all be- 
ing < 0.3 standard deviations from zero. The means for the 
subsamples are also typically < 0.02 mag from zero, with no 
clear bias in any particular subsample. Furthermore, there 
are very few significant outliers in any colour, with only 2- 
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Figure 3. Comparison of synthetic colours derived from our spectra to those measured from light curves at the same epoch. Only spectra 
of objects where there is no obvious galaxy contamination at the position of the SN (as determined from late-time imaging) are included. 
Clockwise from the upper-left panel, we present the B — V,V — R, B — I, and R — I colours. The size of each circle represents the 
size of the photometric uncertainty, with larger circles representing smaller uncertainty. The solid and dotted lines in each panel are the 
one-to-one correspondence and the one-to-one correspondence ±0.1 mag, respectively. 



5 per cent of the spectra (depending on the colour) > 2a 
away from zero. 

In summary, our relative spectrophotometry is excel- 
lent. In particular, objects with little galaxy contamination 
or those where we are able to correct for galaxy contamina- 
tion have extremely good relative spectrophotometry. This 
is achieved simply through our reduction methods and the 
relatively simple host-galaxy contamination correction out- 
lined below; there is no spectral warping of any kind to 
achieve these results. 



3.2.2 Absolute Spectrophotometry 

As mentioned above, there are many achromatic effects 
which can affect our absolute spectrophotometry. We can 
correct for these effects if we have concurrent photometry. 
For these cases, we determined the synthetic photometry of 
our spectra and applied a multiplicative factor to scale the 
synthetic photometry to match our true photometry. This 
scaling is a byproduct of correcting for host-galaxy contam- 
ination, as described in Section 13.31 



3.3 Host-Galaxy Contamination 

SNe generally do not exist in isolation. The vast majority 
occur within galaxies, sometimes close to or on top of com- 
plex regions such as spiral arms or H II regions. With pho- 
tometry, one can correct for this by obtaining a template 
image after the SN has faded (or in some cases, before the 
star explodes), and subtracting the template from the im- 
age with the SN, leaving only the SN (e.g., Filippenko & 
Sargent||1986 1. Although this approach is also feasible with 
spectroscopy (obtaining a spectrum at the position of the 
SN after it has faded), it is not practical. Spectroscopy time 
is typically more valued, and reproducing the exact condi- 
tions at the time of the original SN observation is difficult. 
We do, however, have methods for reducing the galaxy con- 
tamination in a SN spectrum. 

The first method is local background subtraction, as 
described in Section |3.1| Briefly, while extracting the SN 
spectrum, we model the underlying background by interpo- 
lating between background regions on either side of the SN. 
If the background is relatively smooth and monotonic be- 
tween the background regions, this method works very well. 
However, if the SN is near the nucleus of a galaxy or on a 
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Figure 4. Comparison of synthetic colours derived from our spectra to those measured from hght curves at the same epoch. Clockwise 
from the upper-left panel, we present the B — V,V — R, B — I, and R — I colours. The green triangles, blue squares, red diamonds, 
grey circles, and black circles represent (respectively) spectra not observed at the parallactic angle and at high airmass; spectra with 
t > 20 d observed at the parallactic angle, but lacking host-galaxy images to perform host-galaxy subtraction; spectra with t < 20 d 
observed at the parallactic angle, but lacking host-galaxy images to perform host-galaxy subtraction; spectra observed at the parallactic 
angle, with host-galaxy images, but no host-galaxy subtraction is required; and spectra observed at the parallactic angle that have been 
corrected for host-galaxy contamination. The size of each symbol represents the size of the photometric uncertainty, with larger symbols 
representing smaller uncertainty. Representative sizes (0.1 and 0.2 mag) are shown with error bars in the upper- left corner of each 
panel. The short-dashed, dotted, and long-dashed lines in each panel are the one-to-one correspondence, the one-to-one correspondence 
±0.15 mag, and the one-to-one correspondence ±0.3 mag, respectively. 



spiral arm or other bright feature, this method can underes- 
timate the background, leaving galaxy contamination in the 
SN spectrum. 



We have derived a method for removing the residual 
galaxy contamination from our SN spectra. This approach, 
which we call "colour matching," requires both SN photom- 
etry at the time the spectrum was obtained, and template 
colours for the host galaxy at the position of the SN. We 
use the host-galaxy colours to determine the spectral en- 
ergy distribution (SED) of the host galaxy at the position 
of the SN. We then subtract the host-galaxy SED from the 
SN spectrum, scaled so that the synthetic photometry from 
the galaxy-corrected SN spectrum matches the SN photom- 
etry. This method was first presented by |Foley et al.| ( [2010a[ ); 
we discuss it in detail below. 



3.3.1 Determining the Host-Galaxy SED 

The parameter space of galaxy SEDs is well known and well 
behaved, allowing one to reliably reconstruct galaxy SEDs 
with broad-band photometry. Adopting the approach de- 



scribed by Blanton et al. (20031, but updated by Blanton 



& Rowels (20071 to include UV wavelengths, and imple- 



mented in the IDL software package kcorrect . v4_l_4, we 
have used our BVRI photometry of the host galaxy at the 
position of the SN and the redshifts presented in Table [T] to 
reconstruct the galaxy SED at the position of the SN. We 
perform aperture photometry on galaxy templates obtained 
as part of the Lick Observatory SN Search (LOSS) follow- 
up photometry effort ( jGaneshalingam et al.|[2010[ ) using a 
3 pixel (2.4", similar to our Kast slit size and the typical 
seeing at Lick Observatory) aperture and taking the median 
pixel value of the image to represent the sky background. 
Using a 3 pixel aperture for all of our galaxy templates will 
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Table 3. Relative Spectrophotometric Accuracy for the BSNIP Sample 
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50 


S/N < 20 


0.104 


0.093 


0.088 


0.195 


18 


20 ^ S/N < 30 


0.077 


0.026 


0.074 


0.152 


16 


30 ^ S/N < 40 


0.065 


0.055 


0.061 


0.099 


27 


40 ^ S/N < 50 


0.073 
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S/N > 50 
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0.036 
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Reduced by T. Matheson 


0.108 


0.065 


0.115 


0.211 


6 


Reduced by R. Cliornock 


0.065 


0.050 


0.062 


0.101 


34 


Reduced by R. Foley 


0.058 


0.041 


0.040 


0.107 


43 


Reduced by J. Silverman 


0.085 


0.054 


0.081 


0.144 


51 


Reduced by T. Steele 


0.071 


0.022 


0.034 
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represent different physical sizes depending on the distance 
to the galaxy. An aperture significantly different from that 
of the slit combined with the seeing could incorporate flux 
from stellar populations that do not represent the SED of 
the galaxy at the position of SN. As a check on how aperture 
size affects measured galaxy colour, we also used a 4 pixel 
aperture flxed at the SN position. We find excellent agree- 
ment between the colours derived using a 3 pixel aperture 
with a mean difference 5C 0.02 mag. For the typical galaxy 
with z < 0.5 (which includes all redshifts presented here), 
the SEDs are recovered to < 0.02 mag in all filters ( [Blantori] 
et al.|2003||Blanton fc Roweis||2007[ ) . 



3.3.2 Colour Matching 

3.3.2.1 Motivation One approach to subtract galaxy 
contamination from a SN is to extract the SN without any 
local background subtraction, creating a spectrum that con- 
sists of all light at the position of the SN (including galaxy 
light) at the time of the spectrum. If one also has photometry 
at that epoch, one can, in principle, scale a galaxy SED to 
match the galaxy photometry, scale the spectrum to match 
the addition of the SN and galaxy photometry, and sub- 
tract the latter from the former to obtain a SN spectrum. 
The main drawbacks of this method are that (1) one must 
know the proper point-spread function (PSF) of the SN and 
galaxy when the spectrum was obtained, and (2) if there is a 
significant amount of galaxy contamination and the galaxy 
SED is incorrect, significant errors will be introduced. 

When extracting our spectra, we attempt to remove 
as much galaxy contamination as possible. This approach 
has the benefit of reducing the galaxy contamination in the 



SN spectrum without introducing potential errors associ- 
ated with an imprecise photometrically reconstructed galaxy 
SED. Also, considering the lack of precise observing infor- 
mation for many of our spectra (which date back over two 
decades), it would be difficult to estimate the correct PSF 
to determine the exact galaxy flux (both SED and amount) 
entering our slit for a given observation. 

Since the galaxy colours from photometry (which are 
easier to measure than the absolute flux entering our slit) 
determine the galaxy SED, if our spectrophotometry is well- 
calibrated then simply subtracting the galaxy SED until the 
colours of the spectrum match those of the SN photometry 
will result in a SN-only spectrum. 

We can demonstrate this mathematically. In general, an 
observed SN spectrum is defined by 



/spec = A (/sN + B,fg. 



(1) 



where /spec, /sn, and /gai are the vectors of fluxes in the ob- 
served spectrum, SN-only spectrum, and galaxy spectrum, 
respectively, and A and B are normalisation factors. One 
can think of A as normalizing the spectrum in an absolute 
sense to account for slit losses, clouds, and other achromatic 
effects. The parameter B controls the amount of galaxy con- 
tamination, where _B = if there is no galaxy contamination 
and we impose B 0. In principle, B could be negative in 
order to correct for oversubtraction of galaxy light, but our 
testing indicates that allowing B to have negative values 
produces too much overfitting of the spectra. 

From our image templates, we have Pgai, the broad-band 
photometry (in flux units) for the host galaxy at the position 
of the SN. Using MLCS2k2 (|Jha et al.P2"007) template light 



curves or spline interpolations (see Section 3.2.11, we are 
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able to interpolate our SN photometry (independently in 
each band) to determine psN, the broad-band photometry 
(in flux units) for the SN at the time the spectrum was 
obtained. 

We can define the function which translates spectra to 
synthetic broad-band photometry as P, where -P(/sn) = Psn 
and P(/gai) = Pgai- This function is equivalent to convolving 
a spectrum with a filter function. Note that we impose the 
first relationship, while the second relationship is required 
by our method of determining /gai. 

From our spectrum, we are able to determine Pspcc ~ 
P{f spec), the broad-band synthetic photometry (in flux 
units) of the spectrum, which includes both SN and galaxy 
light. These vectors then obey the equation 



-P(/spcc) = A (psN + -Bpgal) ■ 



(2) 



For Equation [2] to be valid, we make two assumptions. 
The first assumption, which is already noted above, is that 
our spectra have accurate relative spectrophotometry. The 
second assumption is that B, the relative fraction of the 
galaxy and SN light, does not vary strongly with wavelength. 
From Section |3.2.1| we have shown that the relative spec- 
trophotometry of our spectra is accurate to ~0. 05-0.1 mag 
across large wavelength regions, comparable to the uncer- 
tainties of our photometry (after interpolating to a given 
date) . 

Solving for /sn in Equation [l] we have 

A /sp, 



/sN 



-B/gal. 



(3) 



With a spectrum spanning at least two bands also covered 
by SN and galaxy photometry, one can solve for A and B 
from Equation [2] With galaxy photometry, the galaxy SED 
(./gal) can be properly reconstructed. It is then simple to 
determine the uncontaminated SN spectrum (/sn) from the 
galaxy-contaminated, observed spectrum (/spec)- We note 
that if B = 0, then Equation |3] simplifies to merely scaling 
the spectrum to match the photometry in an absolute sense. 

3.3.2.2 Testing To test this method, we have performed 
Monte Carlo simulations on six different spectra with in- 
creasing galaxy contamination and appropriate photometric 
errors. Three of the spectra are linear (in fx) and have neg- 
ative, zero, and positive slopes (corresponding to blue, fiat, 
and red spectra). The other three spectra are SN 2005cf 
at maximum brightness, ~ 1 month after maximum, and 
~1 yr after maximum. To each of these spectra we added 
5 galaxy templates, those used by the Sloan Digital Sky 
Survey (SDSS) to perform redshift cross-correlations, span- 
ning early to late galaxy types|^We measured the synthetic 
photometry of the spectra and galaxy templates, and for 
each iteration we varied the photometric data randomly us- 
ing a normal distribution with width corresponding to the 
median error in each band for SNe and galaxies, respec- 
tively. We then performed the colour-matching technique 
for the galaxy-contaminated spectra with the Monte-Carlo- 
based photometry. 

Our recovered SN spectra were compared to our input 
spectra, and the differences between the standard deviation 
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Figure 5. Differences between the median of the standard devia- 
tion of the residuals of contaminated spectra and recovered spec- 
tra after colour matching for different input spectra and galaxy 
templates with varying amounts of galaxy contamination. The 
top through bottom panels correspond to input spectra of a blue 
linear spectrum, a flat linear spectrum, a red linear spectrum, 
SN 2005cf at maximum brightness, SN 2005cf ~ 1 month after 
maximum, and SN 2005cf ^ 1 yr after maximum, respectively. 
Positive values imply that our colour-matching technique yields 
spectra that are closer to the input spectra than the contaminated 
spectra are. The horizontal dotted line in each panel represents 
where the residuals of the recovered spectra are equal to those of 
the contaminated spectra. The blue and red lines correspond to 
the latest and earliest galaxy templates, respectively. The black 
lines correspond to the average over 5 galaxy templates. 



of the residuals of the contaminated spectra and the recov- 
ered spectra are presented in FigurejS] We see that the resid- 
uals for the recovered spectra are significantly lower (i.e., 
the recovered spectra are better at reproducing the input 
spectra) than the contaminated spectra for galaxy contami- 
nations < 70 per cent. The improvement does depend some- 
what on the colour of the SN spectrum and the colour of the 
galaxy template, but the differences are relatively small. At 
higher levels of galaxy contamination, the gains are minimal 
in this metric, but examining the spectra, it is obvious that 
this technique yields impressive results even with very large 
amounts of galaxy contamination. 

In Figure |6j we present our maximum-light and nebular- 
phase spectra of SN 2005cf with varying amounts of galaxy 
contamination. At 70 per cent galaxy contamination, where 
the residuals are not very large, we see that the overall shape 
of the spectra and spectral lines are well recovered. Even at 
90 per cent galaxy contamination, where the contaminated 
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Figure 6. SN spectra used for testing the colour-matching method. The left and right columns correspond to SN 2005cf at maximum 
brightness and ^1 yr after maximum, respectively. Each row represents different amounts of galaxy contamination, from 10 per cent 
(top) to 90 per cent (bottom); each panel has the percent contamination labelled. Each panel shows the input spectra (black), average 
galaxy-contaminated spectra (from several Monte Carlo realisations, red), and average recovered spectra (from several Monte Carlo 
realisations, blue). The maximum-light and nebular-phase spectra have been contaminated with early-type and late- type galaxy templates, 
respectively. The spectra are scaled to have the same median value. Comparing the black (input) spectra to the blue (recovered) spectra 
gives an indication of how well the colour-matching method works. In the top panels, it is difficult to see the input spectra because of 
how closely the recovered spectra match the input spectra. However, even at low levels of contamination, residuals from narrow emission 
lines in the galaxy spectra are seen in the recovered SN spectra. 



spectra appear to be simply galaxy spectra, the method is 
able to recover the overall shape of the SN spectrum. 

The recovered spectra difter most at the ends of the 
spectral coverage. This is due to the galaxy SED reconstruc- 
tion being unconstrained beyond these wavelengths. If we 
extended our galaxy photometry beyond BVRI, this would 
improve. The emission lines of the reconstructed galaxy 
spectra generally have the incorrect strength. This is dif- 



ficult to model with broad-band photometry, and these re- 
gions of the spectra should be ignored. The majority of the 
differences between the input and recovered spectra are the 
result of incorrect galaxy SED reconstruction from errors 
in the galaxy photometry. Improving the galaxy photome- 
try or increasing the number of bands of galaxy photometry 
would improve the reconstruction of the galaxy SED. As 
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the galaxy contamination increases, the errors in the recon- 
structed galaxy SED are amphfied. 

3.3.2.3 Implementation We have applied this tech- 
nique to all SN spectra that have (1) BVRI photometry 
within 5 d of when the spectrum was taken, and (2) a wave- 
length range which spans at least two observed bands. Spec- 
tra which cover only a single observed band are scaled to 
match the photometry at the time of the spectrum. 

The procedure used to subtract galaxy light from an 
observed spectrum is as follows. Using kcorrect, the galaxy 
SED is reconstructed from the broad-band galaxy photome- 
try at the position of the SN. Synthetic photometry is mea- 
sured from the observed spectrum. The SN photometry at 
the time of the spectrum is measured from the light curves 
as described in Section |3.2.1| Using Equation |2] above, the 
factors A and B are determined using a minimisation 
technique. Using Equation [3] the reconstructed galaxy SED 
is subtracted from the observed spectrum to produce the 
corrected SN spectrum. 



4 DATA MANAGEMENT AND STORAGE 

When preparing to present a dataset as large as ours, we 
required some sort of internal organised storage and retrieval 
method. In addition, the overall utility of our dataset will 
be greatly increased by having a user-friendly interface to 
access the data. 

In addition to the final data products, all other informa- 
tion regarding both our photometric and spectroscopic sam- 
ples is stored in our SN Database (SNDB). The SNDB holds 
information about individual SNe (such as host-galaxy in- 
formation, type, discovery information, etc.), much of which 
comes from external, online resources QThe SNDB also con- 
tains information regarding individual spectra (such as ob- 
serving conditions, instrument, resolution, etc.), and indi- 
vidual light curves (number of points, photometric accuracy, 
derived light-curve parameters from various fitting routines, 
etc.). A complete list of all fields stored in the SNDB can be 
found in Table H 

The SNDB contains our entire previously published 
spectral dataset (both SNe la and core-collapse SNe) as well 
as all of the data presented here. It also contains photom- 
etry and light-curve information which has been previously 
published, in addition to photometric data which have been 
compiled and refit by Ganeshalingam et al. (in preparation). 

The SNDB uses the popular open-source software stack 
known as LAMP: the Linux operating system, the Apache 
web server, the MySQL relational database management 
system, and the PHP server-side scripting language. We 
have also implemented instances of the PHP helper classes 
taij^ and JpGraplj^ as well as the JavaScript libraries 

^ For example, lAU Central Bureau for Astronomical Tele- 
grams (|http : //www, cbat ■ eps .harvard, edu/lndex. html] 

fhttp: //www, cbat . eps . harvard. edu/lists/Supernovae . html[l, 

NASA/IPAC Extragalactic Database (NED, |http:| 
[//nedwww. ipac . caltech. edu/ ), and Rochester Academy of 
Sciences Bright Supernova List (http://www.supernovae.net/). 
^ v2.2. Josh Barger (joshbOnpt . com) 

^ vl.27.1, Aditus Consulting ( http : //wwu . aditus . nu/ j pgraph/^ 



SortTabl^and overLI^to improve the functionality and 
user friendliness of the SNDB; we are grateful to the authors 
of these packages. The database is stored on machines at UC 
Berkeley and multiple backups exist at other locations. 

The primary way of accessing the SNDB is via the SN 
Database Public Home Pag^ From here, users can down- 
load precompiled datasets and access our public search page, 
where they can define various input search criteria and query 
the SNDB. All SNe, spectra, and light curves that match 
the search criteria are returned as an HTML table. The re- 
turned information are also written to a WS^^ table which 
is linked from the Search Results Page. If spectra or pho- 
tometry points are returned, users are given the option to 
download the actual data or plot the spectra or photometry 
directly in their web browser. If MLCS2k2 light-curve fits 
are returned, users are given the option to download a file 
containing the fits and probability distributions for each of 
the light-curve parameters. 

As PI Filippenko's group at UC Berkeley publishes 
more spectral and photometric data of SNe of all types in 
the future, the SNDB will continually be updated with these 
newly released data. We hope that the SNDB and its free, 
online access will quickly become an invaluable tool to the 
SN community for the foreseeable future. 



5 CLASSIFICATION 



Optical spectra are often used to classify SNe (e.g., |Filip-| 
penko|1997 Turatto|2003 1 into four basic types. SN II spec- 
tra are identified by strong hydrogen lines which are absent 
in SN I spectra. SN la spectra are characterised by the pres- 
ence of a strong Si II A6355 line typically observed in ab- 
sorption near 6150 A. SN lb spectra lack this Si II feature 
but do contain strong helium lines, and finally SN Ic spec- 
tra lack strong helium lines and have a Si II A6355 line that 
is significantly weaker than those found in SNe la. We per- 
formed automated spectral classification of our full spectral 
datase1p°] using th e SuperNova IDentification code (SNID; 
Blondin & Tonry 20071. Details of our classification algo- 



rithm are presented below. 



5.1 SNID Spectral Templates 

SNID classifies SN spectra by cross-correlating an input 
spectrum with a large database of observed SN spectra 
(known as "templates") which have been deredshifted to 
the rest frame. In order to improve the accuracy of SNID 
classifications, we decided to create our own set of SNID 
spectral templates based on a combination of the default 
SNID templates and our own spectral dataset. 



^ v2, Stuart Langridge jhttp : //www ■ kryogenix ■ org/code/| 
[browser /sorttable/| 



^ v4.21, Erik Bosrup i http: //www. bosrup. com/web/overlib/ I 



http://hercules.b6rkeley.edu/database/lndex_publlc. ] 



Ihtml 



^*~Our full dataset consists of (1) previously published SN spectra 
of all types, (2) SN la spectra which are published here for the 
first time, and (3) some unpublished SN spectra of non-la types 
which will be published in the future. 
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Table 4. SN Database (SNDB) Fields 



SN Name 
Right Ascension 
Declination 
Type 

Type Reference 



SNID-Determined Subtype^ 
Discovery Date 
Discoverer 
Discovery Reference 
Galactic Reddening 



Host-Galaxy Name 
Host-Galaxy Type 
Host-Galaxy Redshift (and error) 
SN Redshift (and error) 
Other Notes 



Photometry and Light-Curve Information 



# of Total Photometry Points 

# of _B-band Photometry Points 

# of V-hanA Photometry Points 

# of _R-band Photometry Points 

# of /-band Photometry Points 

# of Unfiltered Photometry Points 
Photometry Data 

Light Curve Reference(s) 



Ami5(iJ) (and error) 

Maximum S-band Magnitude (and error) 

{B — V) g (and error) 

SALT/2 Distance Modulus (and error)'' 

SALT/2 Light-Curve Stretch (and error)<= 

SALT/2 c (and error)<= 

SALT/2 niB (and error)'' 

SALT/2 xVdof' 



Julian Date of _B-band Maximum (and error) 

Plots of MLCS2k2 Fitsb 

MLCS2k2 Distance Modulus (and error)"^ 

MLCS2k2 Av (and error)*^ 

MLCS2k2 Rv (and error)'' 

MLCS2k2 A (and error)'' 

MLCS2k2 my (and rror)'' 

MLCS2k2 x^/dof'' 



Spectral Information 



# of Spectra of a Given SN 
UT Date of Spectrum 
Filename 

Wavelength Range (A) 

Airmass 

Seeing (arcsec) 

Spectral Resolution(s) (A) 

S/N 



Exposure Time (s) 
Position Angle (deg) 
Parallactic Angle (deg) 
(Observer-Frame) Age 
Observer(s) 
Reducer 
Instrument 
Flux Correction'' 



SNID Redshift (and error)'' 
SNID Type and Subtype'' 

SNID-Determined (Rest-Frame) Age (and error)'' 
SNID rlap" 

SNID Best Matching Template'' 
Spectral Reference(s) 
Flux Standard Star(s) 



"See Section 



''See Jha et a. 



f or mo re information about SNID and its parameters. 

for more information about MLCS2k2 and its parameters. 

and Guy et al. I 2007 1 for more information about SALT and SALT2 and their parameters. 
''See Section 3.3 for more information about our flux corrections. 



See [Guy et al.| 



an| 2007| 



5.1.1 New SNID Subtypes 

We began by downloading SNID vS.cP'] which includes a de- 
fault set of templates consisting of nearby (z < 0.1) SNe of 
all types (la, lb, Ic, II), as well as "NotSN" types, which in- 
cludes galaxies, active galactic nuclei (AGNs), luminous blue 
variables (LBVs), and M stars (see Blondin fc Tonry 2007, 
for the complete default SNID template set). SNID further 
divides each basic SN type into the following subtypes: la- 
norm, Ia-91T, Ia-91bg, la-csm, la-pec, Ib-norm, Ib-pec, lib, 
Ic-norm, Ic-pec, Ic-broad, IIP, Il-pec, Iln, and IIL. 

"Norm" and "pec" identify spectroscopically "normal" 
and "peculiar" SNe of their respective types. Detailed de- 
scriptions of the other subtypes can be found in [Blondin fc| 
|Tonry| ^007\ and [Foley eFaT] (|2009b). In addition to these 
default subtypes, we have added two new SN la subtypes: 
Ia-99aa and Ia-02cx. 

"Ia-99aa" SNe have spectra that resemble those of 
SN 1999aa-Uke obje cts \U et al.[[200T| [Strolger et "ari[2002[ 
Garavini et al.|2004 |. Before maximum brightness, 99aa-like 



SNe contain a Si II A6355 absorption line that is stronger 
than those seen in 91T-like objects, but weaker than those 
of "normal" SNe la. 99aa-like objects also exhibit prominent 
Fe II and Fe HI features at early epochs, similar to the 91T- 



[http: //marww. In2p3 .frZ-blondln/software/snid/ index . [ 
[Etil[ 



like SNe. Moreover, 99aa-like SNe have strong Ca II H&K 
absorption, as do normal SNe la, but quite in contrast with 
91T-like SNe which lack this feature; this is the main spec- 
troscopic difference between 91T-like and 99aa-like objects. 
A comparison of early-time spectra of a 99aa-like SN, a 91T- 
like SN, and a "normal" SN la is shown in Figure [7| SN 
99aa-like events were previously included in the Ia-91T sub- 
type in the default set of SNID templates, but we feel that 
they perhaps represent a spectroscopically distinct subclass 
and therefore deserve their own subtype in SNID. 

Similarly, "Ia-02cx" SNe have spectra that resemble 
those of SN 2002cx-like objec ts (e.g., [Li et al.[[2003[ [Jha[ 
et al.[[2006[ [Foley et al.|[2009a[ ). These SNe were previously 



included in the la-pec subtype in the default set of SNID 
templates, but again, we believe that they represent their 
own subclass of events and should have their own subtype 
in SNID. Note that for our purposes, the "la-pec" category 
refers mainly to SN 2000cx-like objects ( |Li et al.|200T| ). 

We have made a few further changes to the classifica- 
tion scheme of SNID. Namely, Type lib SNe (whose spectra 
evolve from a Type II SN to a Type lb SN, as in SNe 1987K 



and 1993J; see 


Filippenkol |1988| [Filippenko et al.||1993| 


Matheson et al. 


2000 


1 are included only in the "lb" SNID 



type (as opposed to being included in the "11" SNID type 
as well). We have also added two subtypes to the "NotSN" 
SNID type: quasi-stellar objects (QSOs) and carbon stars. 
Spectra of these objects were obtained from the SDSS Data 
Release 6 Spectral Cross-Correlation Templates. 
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4000 5000 6000 7000 ^ 8000 9000 
Rest Wavelength (A) 

Figure 7. Spectra of SN 1991T, SN 1999aa, and the "normal" Type la SN 2006ax at 9.1, 10.6, and 10.1 d before maximum brightness, 
respectively. All spectra in the figure (as well as all spectra plotted in this work) have had their host-galaxy recession velocities removed 
and have been corrected for MW reddening according to the values presented in Tableland assuming that the extinction follows the 
[Cardelli et al.| ( [l989[ l extinction law modified by |0'Donnell| ( [l994^ . Notice how the Si II A6355 absorption feature (marked by the dashed 
line) increases in strength from SN 1991T to SN 1999aa to SN 2006ax. Also note the lack of Ca II H&K absorption (marked by the 
dotted line) in SN 1991T, which is the major difference between 91T-like and 99aa-like SNe. 



5.1.2 New SNID Templates 

We began constructing our new set of spectral templates for 
SNID by performing a literature search for the best-studied 
and most canonical SNe in each subtype (with emphasis 
on the various SN la subtypes). Of the objects we deemed 
to be "the best" examples of their respective subtypes, 30 
are already included in the default set of SNID v5.0 tem- 
plates HBlondin fc Tonry|2007| ) and 30 are in our full spectral 
dataset, with 23 found in both sets. These objects make up 
version 1.0 (vl.O) of our new spectral template set. Tablejs] 
contains information about a subset of the objects included 
in vl.O, as well as the rest of our final template set (the full 
table is available online - see Supporting Information). We 
present a summary of the number of each (sub)type of SN 
included in the final template set in Table [6] Figure |8] shows 
a histogram of the ages of our SN la template spectra. 
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Table 5: SNID v7.0 Spectral Templates 



SN Name 


Subtype 


Version'' 


Age(s)b 


SN 1988Z 


Iln 


1 




SN 1989B<= 


la-norm 


1 


-6, -1,4,6,8, 10,[12:14], [16:25] ,31,37,49,50(1) 


SN 1990H 


IIP 


3 




SN 1990N': 


la-norm 


1 


-13,-6,3,5,15,18,39,50(8) 


SN 1990Q 


IIP 


3 




SN 1991C 


Iln 


2 




SN 1991T'^ 


Ia-91T 


1 


-12,[-10:-5], 0,7,11, 16, 19,25, [46:47], 50(5) 


SN 1991ao 


IIP 


3 




SN 1991av 


Iln 


2 




SN 1991bg'* 


Ia-91bg 


1 


[0:3], [15:16], [19:20], [26:27],30, [33:34], [46:48], 50(10) 


SN 1992A'= 


la-norm 


1 


-5,[-l:0],[2:3],[6:7],9,12,[16:17],24,28 


SN 1992H 


IIP 


3 




SN 1992ad 


IIP 


3 




SN 1993E 


IIP 


3 




SN 1993G 


IIP 


3 




SN 1993J<* 


lib 


1 


[-18:-16], -11, [-5:-2],l, [3:7], [11:13], 17,20,24,[29:30],32,[35:39],41, 50(38) 


SN 1993W 


IIP 


3 




SN 1993ad 


IIP 


3 




SN 1994Dd 


la-norm 


1 


[-12:-2],0,[2:3], [10:17], 19,21,24,26,28,30,40,[43:44],46,[48:49],50(9) 


SN 19941^ 


Ic-norm 


1 


-6, [-4:-3], [0:3], [21:24], 26,30,36,38,40,50(1) 


SN 1994S 


la-norm 


3 


2 


SN 1994W 


Iln 


1 




SN 1994Y 


Iln 


2 




SN 1994ae<^ 


la-norm 


1 


[1:4],6,[9:11],30,36,40,50(6) 


SN 1994ak 


Iln 


3 




SN 1995Dd 


la-norm 


1 


4,6,8,[10:11], 14,16,33,38,43,50(3) 


SN 1995E 


la-norm 


3 


-2 


SN 1995G 


Iln 


2 




SN 1995J 


IIP 


3 




SN 1995V 


IIP 


3 




SN 1995X 


IIP 


3 




SN 1995ac'* 


Ia-91T 


2 


-6,24 


SN 1996L<^ 


Iln 


1 


8,34,42,50(4) 


SN 1996ae 


Iln 


2 




SN 1996an 


IIP 


3 




SN 1996CC 


IIP 


3 




SN 1997Y 


la-norm 


3 


2 


SN 1997ab 


Iln 


2 




SN 1997br<i 


Ia-91T 


2 


[-9:-6],-4,[8:9], 12, 17,21,24,38,42,46,49,50(3) 


SN 1997da 


IIP 


3 




SN 1997dd 


lib 


2 




SN 1997efd 


Ic-broad 


1 


-14,[-12:-9],[-6:-4],7,[13:14], 17,22, 24,27,38,40,44,46,48,50(4) 


SN 1997eg 


Iln 


2 




SN 1998A 


IIP 


3 




SN 1998E 


Iln 


2 




SN 1998S'^ 


Iln 


1 


-13,-2, [1:3], [10:11], [13:14], 16,[31:32],[40:41],44,[46:47],50(37) 


SN 1998bud 


la-norm 


1 


[-3:-l],l,[9:14],[28:44],50(9) 


SN 1998bw= 


Ic-broad 


1 


[8:9], [12:14], 16,[18:19],21,24,[26:28],34,37,43,50(5) 


SN 1998dl 


IIP 


3 





Table abridged; the full table is available online — see Supporting Information. 
All spectral templates are solely from our full dataset, unless otherwise noted. 

"Version of new SNID spectral templates when object was added — 1: vl.O; 2: v2.0— v2.5; 3: v3.0— v7.0. 
'^Rest-frame SN age(s), rounded to nearest whole day, in days from _B-band maximum (for SNe la), from V-band 
maximum (for SNe Ib/c), or from the estimated date of explosion (for SNe II). Ages of spectral templates from our 
dataset are calculated from the light curve references in Table [T] ages from the original SNID v5.0 set of spectral 
templates are from Blondin & Tonry (2007). Adjacent ages are listed in square brackets. Spectra whose age exceeds 
-1-50 days are grouped together and the number of such spectra are noted in parentheses. Many core-collapse SNe 
from our full spectral dataset lack age information (though we require SNe la templates to have age information) . 
"^Spectral templates are from the original SNID v5.0 set of templates only. 

'^Spectral templates are from both our full dataset as well as the original SNID v5.0 set of templates. 
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For each object in vl.O, we examined only spectra that 
had a S/N of at least 15 pixel"^, a minimum wavelength of 
less than 4500 A, and a maximum wavelength of greater than 
7000 A. We also required that each SN la have a date of max- 
imum brightness either from published sources or from Wang 
et al. (in preparation) so we can accurately calculate the 
age of each spectrum. Finally, each spectrum was visually 
inspected by multiple coauthors to be sure they truly rep- 
resented their supposed subtype and were relatively free of 
host-galaxy contamination. If a spectrum passed the quan- 
titative criteria and the by-eye inspection, it was cropped to 
3500-10,000 A (to remove edge artefacts on both ends of the 
spectra). If a spectrum did not cover this entire range, 50 A 
on both ends of the spectrum were removed instead. Finally, 
the cropped spectrum was made into a template; the result 
of this process was vl.O of our new SNID templates. 

To increase the number of SNe in our template set, we 
ran SNID (with our vl.O templates) on our entire spectral 
dataset. To determine subtypes, we followed similar classifi- 



cation criteria to those of Blondin &: Tonry ( 2007 1, requiring 
that the SNID rlap valu^^'^jbe at least 10 and the 3 best- 
matching spectra from SNID all be of the same subtype. We 
also ignored any objects that were classified as "la-norm" or 
"IIP" since we wanted to concentrate only on relatively rare 
subtypes at this point and SNID is somewhat biased toward 
classifying objects as subtypes that have a large number of 
templates (such as "la-norm" and "IIP"). The result of this 
process was v2.0 of our new SNID templates. This process 
was repeated iteratively, running SNID with the previously 
created version of our spectral templates, until no more SNe 
passed all of the criteria. It required 5 additional runs to 
reach this convergence, resulting in v2.5 of our new SNID 
templates. 

We continued creating a new set of spectral templates 
by running SNID (with our v2.5 templates) on our full spec- 
tral dataset. The result of this process was v3.0 of our new 
SNID templates. This process was repeated iteratively, run- 
ning SNID with our previously created version of the spec- 
tral templates, until no more SNe were classified and turned 
into templates. We finished with the creation of v7.0, our 
final set of new SNID templates, which we use to classify 
the remainder of our full spectral dataset (see Table [5] for 
information regarding the entire SNID template set). This 
v7.0 contains 1543 spectra of 277 SNe, of which 779 spectra 
and 134 objects are SNe la. Again, we show a histogram of 
the ages of all of our SN la spectral templates (separated by 
SNID template version) in Figure |8] 



5.1.3 Final Verifications 

As a sanity check, we perform final classification verifica- 
tions. We ran all of our spectra of the objects in v2.5 through 
SNID using our v7.0 templates, and then we ran all of our 
spectra of the objects in v7.0 through SNID (again using 
our v7.0 templates). We compare the (sub)type of the best- 
matching template to the actual (sub)type of the object. 



making sure to ignore all templates of the SN currently be- 
ing inspected. SNe are not used in this process if their spec- 
tra constitute > 15 per cent of all spectra in the object's 
subtype. 

Using only objects from v2.5 (v7.0) we find that SNID, 
using v7.0 of our new templates, is able to correctly clas- 
sify ~97 per cent (~99 per cent) of SN la spectra as one 
of the SN la subtypes and non-la SN spectra as one of the 
non-la SN subtypes. We also find that SNID is able to cor- 
rectly classify ~85 per cent (~95 per cent) of SN la spectra 
with ages ^ 15 d past maximum brightness as the correct 
subtype. 

We compared the classification results of our full dataset 
using our v7.0 templates to the results using the default set 
of SNID templates. The average difference between Zgai (the 
actual redshift of the host galaxy) and zsnid (the redshift of 
the SN as determined by SNID) was found to decrease us- 
ing our v7.0 templates. Furthermore, the discrepancies be- 
tween ti^c (the spectral age derived from photometry) and 
tsNiD (the SNID-determined spectral age) improved drasti- 
cally when using our v7.0 templates. Our template set also 
markedly suppressed the SNID-determined age bias seen 
near maximum brightness and -1-30 d (see Ostman et i 
2011 and Section 5.3 for more on this bias). Finally, with 
the significant increase in the number of templates of pecu- 
liar SNe la and our two additional SN la subtypes, SNID 
(using our v7.0 template set) can distinguish between the 
various spectroscopic subtypes much better than when us- 
ing the default templates. The SNID-determined subtype of 
~15 per cent (~26 per cent) of the spectra (SNe la) presented 
here is different when using SNID with our v7.0 templates 
versus using SNID with the default templates. 



5.2 Classification of Spectra 

Using our v7.0 SNID templates, we can attempt to classify 
all of the spectra in our full dataset using criteria similar 



to those of Miknaitis et al. (20071 and Foley et al. ( 2009b I 



To do this we execute a series of SNID runs to separately 
determine the type, subtype, redshift, and age of the input 
spectrum. For all of the SNID runs we ignore all templates 
of the SN currently being inspected so that SNID will not 
match an object to itself, and we truncate all spectra at 
10,000 A to avoid any possible second-order light contami- 
nation. Besides these, we use the default parameters of SNID 
unless specifically noted below. If a spectrum was obtained 
within 10° of the parallactic angle (or was obtained at an air- 
mass < 1.1) and corrected for host-galaxy contamination via 
our colour-matching technique (as described in Section [3.3[ ), 
then we use the galaxy-subtracted spectrum in all of the 
SNID runs. A subset of the results of the classification al- 
gorithm presented below can be found in Table [t] (the full 
results are available online - see Supporting Information). 



In SNID, the rlap value is a measure of the strength of the 
correlation between the best-matching spectrum and the input 
spectrum. 
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Table 6. Summary of SNID v7.0 Spectral Templates 



la (Total) 


134 


lb (Total) 


16 


Ic (Total) 


14 


II (Total) 


113 


NotSN (Total) 


29 


la-norm 


101 


Ib-norm 


6 


Ic-norm 


10 


IIP 


76 


AGN 


1 


Ia-91T 


3 


Ib-pec 





Ic-pec 


1 


Il-pec 


3 


Gal 


11 


Ia-91bg 


16 


lib 


10 


Ic-broad 


3 


Iln 


34 


LBV 


3 


la-csm 


2 










IIL 





M-star 


7 


la-pec 


1 














QSO 


4 


Ia-99aa 


7 














C-star 


3 



Ia-02cx 4 




Rest-Frame Phase (days) 



Figure 8. A histogram of the ages of our SN la template spectra separated by SNID template version. 



22 Silverman, et al. 



Table 7: SNID Classification Information 



SN Name Type Subtype zsnid'' *SNID rlap Best Match'^ 



SN 


1989A 


la 


la-norm 


0.0087 


(0 


0014) 


83.3 (7.4) 


13 


1 


sn99dk 


(la-norm) 





0088 


(0 


0044) 


71 


83 


(6 


95) 


6 


SN 


1989B 


la 


la-norm 


0.0055 


(0 


0039) 


7.5 (1.4) 


17 


4 


sn05ki 


(la-norm) 





0043 


(0 


0031) 


7 


96 


(1 


49) 


14 


SN 


1989B 


la 


la-norm 


0.0046 


(0 


0041) 


23.6 (-1.0) 


10 


9 


sn99dk 


(la-norm) 





0042 


(0 


0048) 


23 


62 


(0 


68) 


2 


SN 


1989B 


la 


la-norm 


0.0032 


(0 


0015) 


87.3 (13.2) 


19 


5 


sn08ec 


(la-norm) 





0033 


(0 


0027) 


71 


57 


(9 


95) 


16 


SN 


1989B 


la 


la-norm 


0.0029 


(0 


0015) 




15 


7 


sn02fk 


(la-norm) 





0035 


(0 


0035) 


148 


74 


(38 


33) 


9 


SN 


1989M 


la 


la-norm 


0.0014 


(0 


0053) 




17 


2 


sn89B 


(la-norm) 





0043 


(0 


0036) 


—6 


30 


(5 


85) 


9 


SN 


1989M 


la 


la-norm 


0.0022 


(0 


0054) 




13 


1 


sn92A 


(la-norm) 





0065 


(0 


0041) 


6 


30 


(4 


96) 


10 


SN 


1989M 


la 


la-norm 


0.0072 


(0 


0018) 




10 


8 


sn90N 


(la-norm) 





0059 


(0 


0053) 


213 


40 


(0 


00) 


1 


SN 


1989M 


la 


la-norm 


0.0061 


(0 


0000) 




8 





sn90N 


(la-norm) 





0061 


(0 


0073) 


332 


50 


(0 


00) 


1 


SN 


1990G 


la 


la-norm 


0.0365 


(0 


0047) 


9.3 (1.1) 


24 


2 


sn98bu 


(la-norm) 





0359 


(0 


0022) 


9 


30 


(3 


25) 


51 


SN 


1990M 


la 


la-norm 


0.0082 


(0 


0000) 


216.9 (-1.0) 


10 


4 


sn98bu 


(la-norm) 





0082 


(0 


0057) 


208 


00 


(6 


29) 


2 


SN 


1990M 


la 


la-norm 


0.0077 


(0 


0019) 




12 


5 


sn02cs 


(la-norm) 





0091 


(0 


0046) 


31 


28 


(0 


00) 


1 


SN 


1990M 


la 




0.0080 


(0 


0042) 




12 


5 


sn91T 


(Ia-91T) 





0101 


(0 


0042) 


24 


80 


(12 


41) 


201 


SN 


1990M 


la 




0.0081 


(0 


0044) 




13 


3 


sn99aa 


(Ia-99aa) 





0060 


(0 


0042) 


59 


43 


(21 


88) 


172 


SN 


1990M 








































SN 


1990O 


la 


la-norm 


0.0291 


(0 


0025) 


15.1 (1.1) 


25 


6 


sn07qe 


(la-norm) 





0318 


(0 


0024) 


16 


00 


(1 


65) 


22 


SN 


1990O 


la 




0.0287 


(0 


0042) 




19 


9 


sn99dq 


(Ia-99aa) 





0310 


(0 


0026) 


24 


09 


(17 


54) 


275 


SN 


1990O 


la 


la-norm 


0.0299 


(0 


0017) 


55.1 (5.2) 


16 


1 


sn02bo 


(la-norm) 





0284 


(0 


0032) 


56 


35 


(7 


54) 


13 


SN 


1990N 


la 


la-norm 


0.0052 


(0 


0054) 


7.1 (2.8) 


22 


5 


sn95D 


(la-norm) 





0030 


(0 


0022) 


8 


10 


(2 


93) 


22 


SN 


1990N 


la 


la-norm 


0.0033 


(0 


0019) 


18.7 (3.7) 


18 





sn05cf 


(la-norm) 





0019 


(0 


0032) 


18 


69 


(3 


20) 


30 


SN 


1990N 


la 


la-norm 


0.0031 


(0 


0008) 


43.1 (4.3) 


16 


6 


sn99dk 


(la-norm) 





0021 


(0 


0034) 


44 


20 


(7 


14) 


41 


SN 


1990N 








































SN 


1990N 


la 


la-norm 


0.0043 


(0 


0013) 




19 





sn94D 


(la-norm) 





0043 


(0 


0031) 


115 


09 


(41 


77) 


10 


SN 


1990R 


la 


la-norm 


0.0175 


(0 


0010) 


41.3 (1.4) 


17 


5 


sn94D 


(la-norm) 





0178 


(0 


0029) 


43 


20 


(4 


29) 


34 


SN 


1990R 


la 




0.0183 


(0 


0042) 




15 


5 


sn99aa 


(Ia-99aa) 





0332 


(0 


0036) 


26 


20 


(19 


30) 


239 


SN 


1990R 


la 


la-norm 


0.0150 


(0 


0013) 


91.1 (2.8) 


18 


2 


sn94D 


(la-norm) 





0172 


(0 


0030) 


74 


23 


(12 


81) 


14 


SN 


1990Y 


la 


la-norm 


0.0410 


(0 


0013) 


16.1 (2.3) 


12 


7 


sn02bg 


(la-norm) 





0397 


(0 


0045) 


-3 


65 


(7 


14) 


8 


SN 


1991B 


la 


la-norm 


0.0096 


(0 


0021) 


40.5 (-1.0) 


11 


3 


sn05de 


(la-norm) 





0089 


(0 


0045) 


40 


49 


(3 


43) 


2 


SN 


1991B 


la 


la-norm 


0.0098 


(0 


0009) 


57.1 (4.2) 


15 


3 


sn04ey 


(la-norm) 





0097 


(0 


0035) 


51 


49 


(6 


27) 


17 


SN 


1991B 


la 




0.0133 


(0 


0045) 




14 


5 


sn08ds 


(Ia-99aa) 





0109 


(0 


0033) 


63 


44 


(28 


87) 


182 


SN 


1991K 








































SN 


1991K 


la 


la-norm 


0.0189 


(0 


0007) 


91.1 (16.8) 


12 





sn07af 


(la-norm) 





0188 


(0 


0043) 


91 


06 


(15 


17) 


4 


SN 


1991M 


la 


la-norm 


0.0055 


(0 


0010) 


22.5 (-1.0) 


10 


6 


snOlbg 


(la-norm) 





0066 


(0 


0051) 


18 


91 


(2 


52) 


2 


SN 


1991M 


la 


la-norm 


0.0066 


(0 


0014) 




9 


1 


sn05de 


(la-norm) 





0062 


(0 


0056) 


40 


49 


(7 


97) 


89 


SN 


1991M 


la 


la-pec 


0.0096 


(0 


0009) 


146.4 (-1.0) 


14 


1 


snOOcx 


(la-pec) 





0091 


(0 


0045) 


146 


31 


(0 


06) 


2 


SN 


1991M 


la 


la-norm 


0.0098 


(0 


0007) 




10 


1 


sn90N 


(la-norm) 





0089 


(0 


0054) 


213 


40 


(0 


00) 


1 


SN 


19910 


la 


Ia-91bg 


0.0365 


(0 


0028) 




12 


4 


sn06em 


(Ia-91bg) 





0391 


(0 


0060) 


20 


95 


(0 


00) 


1 


SN 


1991S 


la 


la-norm 


0.0556 


(0 


0021) 


38.5 (7.1) 


8 


7 


sn94D 


(la-norm) 





0556 


(0 


0051) 


43 


20 


(11 


50) 


124 


SN 


1991T 


la 


Ia-91T 


0.0029 


(0 


0012) 


-5.6 (-1.0) 


8 


6 


sn95ac 


(Ia-91T) 





0053 


(0 


0057) 


-5 


61 


(2 


39) 


2 


SN 


1991T 


la 


Ia-91T 


0.0037 


(0 


0015) 




14 


5 


sn97br 


(Ia-91T) 





0042 


(0 


0037) 


—7 


40 


(0 


00) 


1 


SN 


1991T 


la 


la-norm 


0.0041 


(0 


0015) 




16 


2 


sn94ae 


(la-norm) 





0052 


(0 


0037) 


9 


40 


(0 


00) 


1 


SN 


1991T 


la 




0.0081 


(0 


0042) 




14 


8 


sn98es 


(Ia-99aa) 





0057 


(0 


0038) 


78 


70 


(28 


43) 


143 


SN 


1991T 


la 


la-norm 


0.0061 


(0 


0011) 


87.3 (6.1) 


14 


6 


sn94D 


(la-norm) 





0070 


(0 


0040) 


87 


19 


(13 


58) 


12 


SN 


1991T 








































SN 


1991T 


la 




0.0082 


(0 


0078) 




23 


1 


sn94D 


(la-norm) 





0080 


(0 


0028) 


115 


09 


(91 


53) 


49 


SN 


1991T 








































SN 


1991T 








































SN 


1991am 


la 


la-norm 


0.0604 


(0 


0022) 


16.4 (2.9) 


16 


6 


sn89B 


(la-norm) 





0602 


(0 


0030) 


13 


40 


(3 


79) 


21 


SN 


1991ak 


la 


la-norm 


0.0114 


(0 


0010) 


41.5 (2.0) 


15 


9 


sn94D 


(la-norm) 





0115 


(0 


0031) 


43 


20 


(4 


26) 


31 


SN 


1991ak 


la 


la-norm 


0.0112 


(0 


0016) 


57.5 (1.5) 


19 


3 


sn02cr 


(la-norm) 





0100 


(0 


0023) 


57 


47 


(3 


98) 


13 


SN 


1991ak 


la 




0.0124 


(0 


0042) 




20 


2 


sn08ds 


(Ia-99aa) 





0104 


(0 


0023) 


63 


44 


(26 


22) 


202 


SN 


1991at 


la 


la-norm 


0.0429 


(0 


0021) 


36.4 (5.1) 


12 


2 


sn04dt 


(la-norm) 





0416 


(0 


0036) 


31 


97 


(8 


38) 


12 


SN 


1991as 


la 




0.0135 


(0 


0042) 




17 


4 


snOOcx 


(la-pec) 





0140 


(0 


0034) 


88 


93 


(77 


90) 


78 


SN 


1991ay 


la 


la-norm 


0.0487 


(0 


0016) 


14.5 (1.4) 


17 


8 


sn95D 


(la-norm) 





0471 


(0 


0033) 


16 


10 


(4 


14) 


31 


SN 


1991bd 


la 


la-norm 


0.0144 


(0 


0044) 


23.6 (-1.0) 


11 


6 


sn99dk 


(la-norm) 





0139 


(0 


0049) 


23 


62 


(1 


36) 


2 


SN 


1991bc 


la 


la-norm 


0.0223 


(0 


0026) 


17.7 (2.9) 


24 





sn90N 


(la-norm) 





0219 


(0 


0024) 


17 


70 


(4 


24) 


33 


SN 


1991bc 


la 




0.0233 


(0 


0039) 




14 


9 


sn91T 


(Ia-91T) 





0202 


(0 


0035) 


75 


20 


(22 


67) 


259 



Table abridged; the full table is available online — see Supporting Information. 
The entries in this table match one-to-one with the entries in Table [2] 
"The redshift uncertainty is in parentheses. 

''Phases of spectra are in rest-frame days. The phase uncertainty is in parentheses. Phase uncertainties of imply that only one 
template was a "good" match. 

'^The best matching SNID template in the form: "template SN" (subtype) ZSNID (error) Jsnid (error) N_good_matches 
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5.2.1 SNID Type 



5.3 Verifying Redshifts and Ages from SNID 



If the object's redshift is known a priori (from the host 
galaxy, usually via NED), we force SNID to use this red- 
shift; otherwise we do not use any redshift prior. For the 
first attempt to determine a type, the minimum rlap value 
is set to 10. A spectrum is determined to be of a given type 
if the fraction of "good" correlations that correspond to this 
type exceeds 50 per cent. In addition, we require the best- 
matching SN template to be of this same type. If the spec- 
trum's type cannot be determined using these criteria we 
perform another SNID run, this time using a minimum rlap 
value of 5. This resulted in 1232 of our 1298 spectra receiving 
a SNID type. If the type of the input spectrum is successfully 
determined (using either rlap value), an attempt is made to 
determine its subtype. 

5.2.2 SNID Subtype 

Again, we adopt the host-galaxy redshift when available and 
the minimum rlap used to determine if the subtype is the 
same as was used to successfully determine the type (either 5 
or 10). We also force SNID to only consider templates of the 
previously determined SN type. A spectrum is determined 
to be of a given subtype if the fraction of "good" correla- 
tions that correspond to this subtype exceeds 50 per cent. In 
addition, we require the best-matching SN template to be 
of this same subtype. If the spectrum's subtype cannot be 
determined using these criteria, and a minimum rlap of 10 
was used, we perform another SNID run, this time using a 
minimum rlap value of 5. This resulted in 1098 of our 1298 
spectra receiving a SNID subtype. Regardless of whether 
SNID determines a subtype, a third run is executed to de- 
termine the redshift. 



5.2.3 SNID Redshift 

The SNID redshift is calculated by taking the median of 
all "good" template redshifts and the redshift error is the 
standard deviation of these redshifts. If a subtype has been 
successfully determined, we force SNID to only use tem- 
plates of that subtype; otherwise we only use templates of 
the previously determined type. For this SNID run no a pri- 
ori redshift information is used. 1232 of our 1298 spectra 
received a SNID redshift. If a redshift is successfully deter- 
mined in this run, a fourth run is executed to calculate the 
age of the spectrum. 



5.2.4 SNID Age 

The SNID age is calculated by taking the median of all 
"good" template ages that have an rlap value larger than 
75 per cent of the rlap value of the best-matching template. 
The age error is the standard deviation of these ages. Only 
if a subtype has been successfully determined do we attempt 
to calculate an age. We also require that the age error be 
less than either 4 d or 20 per cent of the SNID-determined 
age (whichever is larger). For this run, we again force SNID 
to adopt the host-galaxy redshift when available. If no host- 
galaxy redshift is known, we use the previously determined 
SNID redshift instead. 849 of our 1298 spectra received a 
SNID age. 



It has been shown previously that SNID-determined red- 
shifts correlate extremely well with actual redshifts of SN 
host galaxies (e.g., Ma t heson et _al.| |2005| [Foley et al.| 



2009b 



Ostman et al.|2011[ ). Figure [g] shows that the SNID- 



determined redshifts (using v7.0 of our templates) agree well 
with the host-galaxy redshifts of our data. The dispersion 
about the one-to-one correspondence is only ~0.002 for the 
1184 spectra for which the redshift is known and which SNID 
determined were SNe la and calculated a redshift. This is as 
good as or better than what has been found previously us- 
ing rnuch smaller samples of higher-redshift SNe ([Matheson 



et al. 



2005 



Foley et al. 



2009b 



Ostman et al. 



20111. How- 



ever, as seen in Figure |9 the majority of the largest outliers 
appear to have SNID redshifts that are lower than the host- 
galaxy redshifts. The normalised median absolute deviation 



(i.e., a measure of the precision of our redshifts; Ilbert et al 



20061, defined as 



1.48 X median 



I^SNID — 2gal| 
1 -I- 2gal 



(4) 



is 0.002. 

One spectrum is a fairly significant outlier (seen in the 
top left of the top panel of Fig. |9|. This spectrum is 360 d 
past maximum brightness (according to the light curve) and 
we only have a small number of SN la templates that are 
this old (only 3 older than 300 d). The relative lack of good 
matches with old SN la spectra, as opposed to much younger 
spectra at much higher redshifts, is most likely the cause of 
the erroneous redshift (and age) from SNID. However, aside 
from this extreme case, the SNID-determined redshifts are 
quite accurate. Only 8 per cent (3 per cent) of the objects 
are more than 2cj (3cr) away from the one-to-one correspon- 
dence. 

The original SNID template spectra have ages which 
have been corrected for the 1/(1 -I- z) time-dilation factor 
that we expect to observe in an expanding universe (e.g., 



Riess et al. 1997 Foley et al. 2005 Blondin et al. 20081 



Thus, SNID templates should have ages in the rest frame 
of the SNel^We compare the SNID-determined ages of our 
SN la spectra to their actual (rest-frame) ages as derived 
from their light curves and redshifts as presented in Table [l] 
the result is shown in Figure [lO] There are 595 total spectra 
that SNID determined were SNe la (and 409 with light- 
curve ages 5C 30 d) which have both SNID-determined ages 
and light-curve ages. The dispersion about the one-to-one 
correspondence for the total sample is ~4.1 d, and the dis- 
persion for the sample with light-curve ages 30 d is ~3.3 d. 



Foley et al. (2009b I calculated a dispersion of 1.8 d for 59 
SN la spectra with light-curve ages between —11 and 19.4 d 
at moderate to high redshift (0.100 < ^ < 0.807). We have 
338 spectra in t his range and c a lculat e a dispersion of ~2.8 d 
for these data. 



Ostman et al. 



(20111 obtained a dispersion 



of 4 d for 127 SN la spectra with light-curve ages between 
-11.7 and 67.9 d at moderate redshift (0.03 z 0.32). 
We have 529 spectra in this range and calculate a dispersion 
of ~3.7 d for these data. 



^•^ When creating our own SNID templates, we transformed our 
SN ages to the rest frame (using the redshift of the SN or the 



host galaxy). 
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Figure 9. Comparison of galaxy redshifts (zgai) and SNID-determined redshifts (zsnid) using our SNID classification sclieme: all 1184 
spectra for which the redshift is known and which SNID determined were SNe la (top) and a zoom-in on the 1114 SN la spectra with 
Zgai 0.05 and ^snid 0.05 (bottom). The median error bar in both directions for the entire sample is shown in the bottom right of 
each plot. The bottom of each plot shows the residuals versus Zgai. In all plots the short-dashed line is the one-to-one correspondence 
and the dotted lines are ^SNID = -^gal i 0.002 (l -|- Zgai) (representing the average precision). See the text for more information on the 
extreme outlier in the top left of the top panel. 
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Figure 10. Comparison of rest-frame light-curve ages (tLc) ™d SNID-determined ages (tsNlo) using our SNID classification scheme; 
all 595 spectra which SNID determined were SNe la and that have both SNID-determined ages and light-curve ages (top), and a zoom-in 
on the 394 SN la spectra with ^! 50 d and tsNiD ^ 30 d (bottom). The median error bar in both directions for the entire sample is 
shown in the top left of each plot. The bottom of both plots shows the residuals versus Jlc- 
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The samples used by Foley et al. ( 2009b I and Ostman 



et al. (2011 1 were at higher redshift than our data, and the 



way in which the y determined spectral ages using SNID was 



Ostman et al. 



(2011 1 simply used the me- 



slightly different, 
dian of all "good" template ages as the SNID-determined 
age. We initially attempted this relatively straightforward 
method for our spectra, but we soon found a significant 
bias in our SNID-determined ages (as compared to ages 
derived from the SN light curves). The bias was caus- 
ing SNID-determined ages to be artificially skewed toward 
about +30 d for spectra which are (according to their pho- 
tometry) in the range of about 23-33 d. We also observed 
a similar, yet weaker, version of this bias for spectra whose 
photometric ages were near maximum brightness as well as 
spectra with light-curve ages near ~100 d (still visible in the 
top panel of Fig. 10 1. The bias near maximum brightncssjs 



seen in the higher-redshift data presented by |Ostman et al.| 
(20111, and we suggest that the bias near ages of ~30 d 



(and perhaps the one near 100 d as well) would have been 
observed in their data if their dataset had contained spectra 
near these epochs. 

We thoroughly investigated these so-called "age attrac- 
tors" and their effect on the dispersions, but unfortunately 
found no simple explanation for the apparent bias. In order 
to characterise and explain the strongest bias (near photo- 
metric ages of one month past maximum brightness), we 
used the default set of SNID templates, instead of our own, 
and again ran all of our spectra through our classification 
routine, only to have the bias appear even stronger than be- 
fore. We examined the age and light-curve shape (as char- 



acterised by the MLCS2k2 A parameter; Jha et al. 20071 



distributions of our v7.0 templates, the default set of SNID 
templates, and the light-curve ages of our entire dataset, 
and none of these showed any deviation at or near -1-30 d 
that might affect SNID's age determinations. For example. 
Figure [8] shows a histogram of the ages of all of our SN la 
spectral templates (separated by SNID template version), 
and there does not appear to be any obvious bias near 
-1-30 d. Furthermore, we investigated how the strength of the 
bias changed with SNID-determined subtype, best-matching 
subtype, and rlap, but saw no strong correlations. 

On the other hand, all reported SNID-determined ages 



from Foley et al. ( 2009b I required that there were at least 



8 good matches (and the SNID-determined age was calcu- 
lated from the median age of those top 8 matches) , that the 
age error (defined as the standard deviation of the top 8 
matches) was < 6 d, and that there was a SNID-determined 
subtype (S. Blondin, 2011, private communication). We also 
attempted to use this method of SNID age determination on 
the BSNIP data but found that an extremely large fraction 
of the spectra were not being assigned a SNID age. 

Various other methods of SNID age determination and 
verification were tested (see Table [s] for a summary), and 
the method that gave the best compromise between low dis- 
persion values and large numbers of spectra being assigned 
a SNID age is the one that we ultimately used (described 
above in Section |5.2.4[ ). Even though this new method of 
determining SNID ages reduces the bias near -1-30 d, the 
bi as near maximum bri ghtness (mentioned above and seen 



Ostman et al. 



of our SNID-derived ages are within la of the one-to-one 
correspondence with light-curve age. 



5.4 Classification of Objects 

After classifying (or attempting to classify) all of the spectra 
in our dataset using the method described in Section [5.2[ we 
use the SNID information for all spectra of a given object 
to determine the SN's (sub)type. 



5.4-1 Classification Accuracy 

To investigate the accuracy of our SN la classification 
scheme, we attempt to find correlations between the accu- 
racy of our SN la subtype determination and the rlap value 
of the best-match template, the SNID-determined age of the 
input spectrum, and the S/N of the input spectrum. We find 
that our accuracy is similarly correlated with both rlap and 
S/N, which is unsurprising since S/N is one of the factors 
which determines the rlap value during a SNID run. Thus, 
we attempt to correlate our classification accuracy with only 
rlap and spectral age, simultaneously. 

To do this, we compare the subtype we determined us- 



ing our classification scheme (Section 5.2 1 for each SN la 
template spectrum in v7.0 to the actual subtype of the tem- 
plate object. If our multiple SNID runs correctly classified a 
given spectrum, we assigned it an "accuracy percentage" 
(P) of 1, and if it was misclassified it received an accu- 
racy classification of 0. We then used the rlap value of the 
best-matching template (rlap) and the SNID-determined 
(rest-frame) age (t, in days) of each spectrum to fit a two- 
dimensional surface of the form 



P = ci + C2 X t + C3 X rlap 



(5) 



20111 is still present. It should also be 



noted that despite the biases that remain, over 60 per cent 



to the SNID classifications of our v7.0 SN la template 
spectra. This function was fit on a grid of phases from 
—20 to 200 d (in steps of 2 d) and rlap values from 5 
to 40 (in steps of 1). Our best-fitting values for the con- 
stants were ci = 0.68 ± 0.04, C2 = -0.0014 ± 0.0002, and 
C3 — 0.018 ± 0.002, and our resulting contours are shown in 
Figure |11| The few spectra which are misclassified but lie 
above the 100 per cent contour are simply ones where even 
though they are misclassified, their age and rlap values for- 
mally imply an "accuracy percentage" of > 100 per cent. 
Using our values for ci, C2, and C3, we can now calculate an 
"accuracy percentage" (P) for any SN la spectrum that has 
a SNID-determined age and rlap (by using Equation [5|. 

As expected, our accuracy increases with increased rlap 
value since rlap is a measure of the strength of the corre- 
lation of the input spectrum with the best-matching SNID 
template. In addition, our accuracy decreases with increased 
age since, as noted earlier, as some subtypes of SNe la evolve 
and age their optical spectra begin to resemble those of "nor- 
mal SNe la." 



5.4-2 Final Object Classification 

For SNe with multiple spectra, we must consider each spec- 
trum's classification when obtaining a final classification for 
a given object. To do this we first determine an object's type 
by counting the total number of spectra of each type for the 
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Table 8. Number of Spectra and Dispersion of Various Methods for Determining SNID Ages 



BSNIP BSNIP BSNIP BSNIP 
(final)'' (alternate)'' (011)'= (F09)<^ 



Oil" 



Fogf 



Age Range (d) 


N 


a (d) 


N 


a (d) 


N 


<x(d) 


N 


a (A) 


N 


<7 (d) 


N 


a (d) 


-20 ^ t ^ 360 


595 


4.1 


869 


5.0 


860 


5.4 


315 


3.2 










-20 ^ t ^ 30 


409 


3.3 


580 


3.9 


572 


4.1 


264 


2.9 










-12 ^ i ^ 68 


529 


3.7 


757 


4.5 


748 


4.7 


314 


3.2 


127 


4 






-11 ^ t ^ 19 


338 


2.8 


451 


3.3 


444 


3.4 


233 


2.4 






59 


1.8 



"Oil" 



Ostman et al. 



''The method described in Section |5.2.4[ 



(20111, "F09" = Foley et al. ( 2009b i 



''Similar to the method described in Section |5. 2. 4| The SNID age is calculated by taking 
the median of all "good" template ages that have an rlap value larger than 75 per cent of 
the rlap value of the b est-matching template . However, no other requirements are imposed. 
''The method used by 



Ostman et al. 



(2011 1 applied to the BSNIP data. This simply uses 
the median of all "goo d" template ages as t he SNID-determined age. 

''The method used by [Foley et al.| ( 2009bl l applied to the BSNIP data. It requires that 
there are at least 8 good matches (and the SNID-determined age is calculated from the 
median age of those top 8 matches), that the age error (defined as the standard deviation 
of the top 8 matches) is < 6 d, and that there is a SNID-determined subtype (S. Blondin, 
2011, private comm unication) . 



"Values reported by [Ostman et al. 



I 2011 using data with 0.03 s£ ^ ^ 0.32. 



Values reported by [I'oley eFaT I pOofE I using data with 0.10 ^ z ^ 0.81 




50 100 150 

Phase (days) 



Figure 11. Accuracy of our SNID-determined SN la subtypes (using our classification scheme described in Section |5.2[ l versus SNID- 
determined (rest-frame) phase (in days) and rlap. The □ are correctly classified spectra from our v7.0 SNID templates; X are misclassified 
spectra from our v7.0 SNID templates. The contours (from bottom to top) represent 60 per cent, 70 per cent, 80 per cent, 90 per cent, 
and 100 per cent accuracy. The few spectra which are misclassified but above the 100 per cent contour are due to the fact that we fit all 
v7.0 template spectra and thus a handful of misclassified objects formally have accuracies above 100 per cent. 
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given object. The object is then classified as the type with 
the highest count. If there is more than one type with the 
same highest count, we compare the spectra of those types 
only and use the type of the spectrum with the largest rlap 
(though we add a "?" to the type to denote our uncertainty 
regarding the classification). 

For each object whose definite type has been determined 
(i.e., no trailing "?") and that is not a SN la, we assign 
a subtype by counting the total number of spectra of each 
subtype for the given object. The SN is then classified as the 
subtype with the highest count. If there is more than one 
subtype with the same highest count, we cannot accurately 
determine the subtype and thus classify the object as simply 
the previously determined type. 

For each object that we have determined is a definite 
SN la (i.e., no trailing "?"), we calculate the "accuracy per- 
centage" (P) for each spectrum of that object using Equa- 
tion|5] We then combine the "accuracy percentages" of all of 
the spectra of the given SN to calculate the maximum prob- 
ability that it is of each subtype. For example, if a spectrum 
has an "accuracy percentage" of P, then that is a measure 
of the probability that the spectrum is of subtype m (where 
subtype m is the subtype determined by our classification 
scheme in Section 5.2 I. This means that 1 — P is the proba- 
bility that this spectrum is any of the other SN la subtypes, 
and it is in fact the maximum probability that the spectrum 
is of subtype n for any n ^ m. Therefore, we can combine 
the "accuracy percentages" for multiple spectra of different 
subtypes to calculate a maximum probability for each sub- 
type. 



Table 9. Summary of Final SN Classifications 



(6) 



where Pm is the maximum probability the given object is 
of subtype m, the first product is over all spectra whose 
subtype (n) is equal to m, the second product is over all 
spectra whose subtype (n) is not equal to m, and km is a 
normalisation constant for subtype m defined as 



(7) 



Y[ Pnx Yi {1 p,o + n (1 - 

n — m n^m n — m 

Once we calculate the maximum probability for each sub- 
type, we classify the given object as the subtype with the 
largest such probability. 

SNID is merely a tool, albeit a useful one for determin- 
ing spectral subtypes of SNe, but ultimately humans classify 
SN spectra. Thus, we visually inspected any objects in our 
dataset that were classified as "la" or "la?" (with no sub- 
type) or as any of the non-SN la subtypes, in addition to any 
object that had spectra that were classified as more than one 
SN la subtype. From these visual inspections we changed 
a handful of the final object classifications from what our 
SNID-based classification algorithm would have yielded. We 
also forced all objects that were v7.0 SNID templates to be 
their actual subtype. These final (sub)type determinations 
can be found in Table [l] we also present a summary of the 
final (sub)type determinations from SNID, as well as our 
adjusted classifications, in Table [9] 

Of all of the v7.0 template objects (which include SNe la 
as well as other SN types), 95 per cent were classified as the 
correct type by SNID and 86 per cent were classified as the 
correct subtype by SNID. We were unable to determine a 
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type (subtype) for 2 per cent (7 per cent) of the v7.0 tem- 
plate objects using our aforementioned SNID-based classifi- 
cation scheme. Of solely the v7.0 template objects that are 
SNe la, 99 per cent received the correct type classification 
from SNID and 92 per cent the correct subtype classifica- 
tion, and we could determine a type and subtype from SNID 
for all of the v7.0 template SNe la. These are more infor- 
mative and more accurate percentages than those from the 
simple sanity checks that were discussed in Section [5.1.3| 

Of the 582 SNe la that are presented in this paper, our 
SNID-based algorithm classifies 96.8 per cent of them as 
SNe la, 1.5 per cent as other SN types, and we are unable to 
classify 1.7 per cent. Over half of the spectra of objects that 
are classified by SNID as non-la types or that are unclassified 
are relatively old and/or noisy, and thus it is not surprising 
that our classification scheme failed on these observations. 
A few spectra of these objects are heavily contaminated by 
host-galaxy light and lacked sufficient photometry for us to 
apply our galaxy-correction algorithm (Section |3.3[ ), and so 
again it is reasonable that these are not correctly classi- 
fied. When it was not clear from our own spectra that these 
objects were SNe la, we relied on previous spectral classifi- 
cations from the literature. In addition, some of the objects 
which were classified as non-la types are well-known, pecu- 
liar SNe la (including a few of our v7.0 SNID templates), 
but we do not have enough templates of other objects of 
their subtype at similar ages for SNID to get good matches 
to our observations. This is most likely responsible for many 
of SNID's "la" (with no definitive subtype) and "la?" classi- 
fications as well, and is partially why we visually reinspected 
these objects. 



5.4-3 Relative Fractions of Subtypes 

Even though our SN la spectral sample is not complete by 
any rigorous definition and suffers from multiple observa- 
tional biases, it is still illuminating to calculate the relative 
fractions of SN la subtypes (as determined by SNID) in our 
dataset. Of the objects for which SNID determines a SN la 
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subtype, 92.1 per cent are normal, 5.4 per cent are 91bg- 
like, 0.6 per cent are 91T-like, 0.4 per cent are 99aa-like, 
and 1.5 per cent are 02cx-like. These fractions are listed in 
Table [To] along with those from other SN la studies. For the 
purposes of comparing these percentages to values in the 



literature, we will follow Li et al. (2011 1 and combine 99aa- 



like objects with 91T-like objects, yielding 1.0 per cent for 
this group. If, as we did when creating our new SNID tem- 
plates in order to get more accurate subtype classifications, 
we require spectra to have ages less than 15 d past maxi- 
mum brightness and S/N > 15 pixel" ^, then we find that 
of these objects 90.7 per cent are normal, 6.5 per cent are 
91bg-like, 1.9 per cent are 91T/99aa-like, and 0.9 per cent 
are 02cx-like. 

In the volume-limited sample of Li et al. (20111, 



70 per cent of SN la were normal, while 15 per cent were 
91bg-like, 9 per cent were 91T/99aa-like (only a lower limit, 
however), and 5 per cent wer e 02cx-like. The sligh tly higher 
redshift {z = 0.17) sample of 



Ostman et al. 



(20111 detected 



only a few probable peculiar SNe la, but they calculate that 
7-32 per cent of local SNe la should be 91T/99aa-like. The 
even higher redshift {z x 0.35) sample of|Foley et al. ( 2009b I 



found that 4-19 per cent were 91T/99aa-like, while they had 
no 91bg-like objects in their dataset However, an analysis 
of the companion photometry to much of the spectroscopic 



sample presented here ( Ganeshalingam et al.|,2010) finds 
(using SNID as described in this work) that 80.6 per cent 
of their SNe la are normal, 10.9 per cent are 91bg-like, 
5.5 per cent are 91T/99aa-like, 2.4 per cent are 02cx-like, 
and 0.6 per cent are la-pec (i.e., SN 2000cx, which SNID 
classified herein as a la-norm since we have no other SNID 
templates of OOcx-like objects). These percentages are closer 
to our fractions of subtypes than the complete sample pre- 



sented by Li et al. (20111. This likely comes from the fact 



that for most of the project's lifetime, BSNIP has had a 
strong focus on spectroscopically monitoring SNe la that 
were being concurrently observed photometrically as part of 
the sample presented by [Ganeshalingam et al.| ( |2010| |. 

The main differences between our SNID-determined 
fractions of SN la subtypes and those found in the complete 



sample of Li et al. ( 2011 1 is that we classify too many objects 



as normal and not enough as 91T/99aa-like and 91bg-like. 
This can possibly be explained by the fact that spectra of 
these peculiar SNe la (especially 91T/99aa-like objects) re- 
semble spectra of normal SNe la within a week or two after 
maximum brightness (introducing the so-called "age bias"; 
Li et al. 20011. Thus, some of our la-norm objects may in 



fact be 91T/99aa-like SNe la, but the spectra in our dataset 
were obtained at too late an epoch to distinguish between 
the two subtypes. Since there are many more normal SNID 
templates than peculiar templates, these objects will ulti- 
mately get classified as normal. Furthermore, it is possible 
that some objects have essentially normal spectra and are 
classified as la-norm by SNID, but have slowly (quickly) de- 
clining light curves and are thus photometrically classified as 



91T/99aa-like (91bg-like) by Li et al. (20111. This interest- 



1* [Foley et af] |2009b| , after correcting for various biases, ex- 
pected only 1—4 91bg-like objects based on the [Li et ah] | |2001| | 
percentages. 



ing possibility will be investigated further in future BSNIP 
papers. 

Many of the non-normal SNe la that we classify in our 
dataset have already been noted as peculiar in the liter- 
ature (either in unrefereed circulars or published papers). 
However, there are still a few peculiar classifications that 
we publish here for the first time. Details regarding these 
interesting individual objects can be found in Section [6.2[ 



6 THE BSNIP SAMPLE 

Our SN la spectral dataset consists of a total of 1298 spectra 
of 582 SNe la observed from 1989 through the end of 2008, 
representing nearly 470 hr of telescope time. 1159 spectra 
of 563 objects are published here for the first time. In Fig- 
ures [12] through [16| we show examples of our most recent, 
fully reduced spectra from each of our commonly used in- 
struments (Kast, LRIS, UV Schmidt, DEIMOS, and ESI, 
respectively). Plots of all of the fully reduced spectra as well 
as (for the objects with multi-band SN and galaxy photome- 



try) galaxy-subtracted spectra (as discussed in Section 3.3 1 
presented in this work are available online - see Support- 
ing Information. Figure [TTj shows an example of a spectral 
sequence, obtained for SN 20071e (spectral sequences for all 
objects in our dataset with more than 7 spectra can be found 
online - see Supporting Information). 



6.1 Sample Characteristics 

If we remove objects where we have no light-curve informa- 
tion, leaving only spectra with phase (and light-curve shape) 
information, the sample is reduced to 914 (770) spectra of 
321 (251) SNe la. If we further restrict ourselves to objects 
with relatively well-sampled, filtered light curves, retain- 
ing only those objects with precise distance measurements, 
our dataset contains 584 spectra of 199 SNe la. Finally, if 
we only count spectra which have reasonable estimates of 
multi-filtered SN magnitudes at the time of the spectrum 
and measurements of the host-galaxy colours at the posi- 
tion of the SN, providing accurate flux-calibrated, galaxy- 



subtracted spectra (see Section 3.3 for more information), 
our sample contains 234 spectra of 95 SNe la. 

This dataset has ~3 times the number of spectra and 
~18 times the number of objects as the sample of [Mathe^ 
[son et aT] ( [2008[ ). Due to the scheduling of their telescope 
time, their dataset consisted mainly of well-sampled spec- 
troscopic time series of a handful of SNe la, averaging 13.5 
spectra per object. By contrast, the BSNIP sample consists 
of ~2.2 spectra per object, showing our emphasis on the to- 
tal number of objects rather than the number of spectra per 
object. The histogram of the number of spectra per object 
for our sample can be seen in Figure [Ts] Thus, our sample 
(emphasizing the total number of objects) and the sample of 
Matheson et al. ( 2008 1 (emphasizing the number of spectra 



per object) are complementary. Furthermore, as mentioned 
in Section[2] our spectra typically cover 3300-10,400 A, com- 
pared to the typical 3700-7400 A wavelength range of the 



spectra from Matheson et al. ( 2008 1 



All of the SNe la in our dataset have z ^ 0.2, and the 
vast majority (89.6 per cent) have z ^ 0.05. The distribu- 
tion of redshifts for all of our SNe la (as well as just those 
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Table 10. Percentages of SN la Subtypes 
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Figure 12. Spectrum of SN 2008hy taken on 2008 Dec. 31 (UT dates are used throughout this Paper) using Kast. This spectrum is 
classified as "la-norm" by SNID. Clipped bad residuals along with host-galaxy emission lines and telluric absorptions (if present) are 
plotted in light grey. 



with phase information) is shown in Figure 19 The aver- 



age redshift of the full sample of objects presented here is 
about 0.0283 and the median uncertainty is 0.00004 (as de- 
termined from the redshift uncertainties reported in NED). 
About 78 per cent of our SNe la have z 0.015 (which is 
approximately the redshift above which peculiar velocities 
can be ignored in cosmological calculations; see, e.g., lAstier] 



et al. 20061. 17 of our objects have unknown host-galaxy 



redshifts, but all of the spectra of these objects received a 



redshift from our SNID classification scheme (Section 5.21. 

As mentioned above, many of our spectra have phase 
information from photometric observations; the distribution 
of phases for all of our spectra is shown in Figure |20[ We 
have 147 spectra of 114 SNe la before maximum brightness. 
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Figure 13. Spectrum of SN 2008ge taken on 2008 Oct. 28 using LRIS. Tliis spectrum is classified as "Ia-02cx" by SNID and "Ia-02cx" 
by [Foley et al.| ( |2010b[ |. Clipped bad residuals along with host-galaxy emission lines and telluric absorptions (if present) are plotted in 
light grey. 



and 245 (107) spectra of 181 (96) objects within 7 d (3 d) of 
maximum brightness. Our sample also contains 34 spectra 
of 20 SNe la older than 180 d past maximum brightness. 
The median uncertainty of our phases is 0.38 d, though the 



practical uncertainty is more like 0.5 d ( Ganeshalingam et al. 
|2010[ ). Figure [21] shows the distribution of the phase of each 
SN la at the time of its first spectrum. 

As mentioned above, many of our SNe have light-curve 
shape information (as characterised here by the MLCS2k2 A 
parameter; Jha et al.|2007 1; the distribution of A values for 
all of these objects is shown in Figure [22] In the histogram 
we also denote each object's final SNID classification (Sec- 
tion 5.4.21. The A values come from previously published 



photometric data which have all been compiled and fit by 
Ganeshalingam et al. (in preparation). The average A for 
our dataset is ~0.12, and the median A and uncertainty are 
about —0.03 and 0.035, respectively. Our SNe la span most 
of the standard range of A values (about —0.4 to 1.6, e.g., 
Hicken et al.|2009 l. 

Figure[23]presents the host-galaxy redshift of our SNe la 
versus their A values. Our dataset covers most of the range 
of A values at the lowest redshifts, but our coverage de- 
creases at higher redshifts. SNe la with large A are the 



fainter, faster-evolving objects (Jha et al. 20071, and thus 



we have fewer of those (relative to smaller A objects) in 
our sample at higher redshifts. Figure [24] presents the phase 
of each spectrum versus their A values. Analogous to Fig- 
ure |23[ our sample spans most of the standard range of A 
values at early times, but the coverage begins to drop off at 



~40 d past maximum brightness. Once again, we lack large- 
A objects at the latest times, while we still have a handful 
of SNe la with small A values at these epochs. This is un- 
surprising since objects with large A values are fainter and 
have faster-evolving light curves, and thus we are not able 
to follow them spectroscopically for as long as their small-A 
brethren. It is interesting to note that there are relatively 
few SNe la with 0.5 sC A ^ 0.7 or 1.1 ^ A sC 1.3; these pos- 
sible anomalies will be investigated further in future BSNIP 
papers. 
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4000 5000 6000 7000 8000 9000 
Rest Wavelength (A) 

Figure 14. Spectrum of SN 1991bg taken 27 d (rest frame) past maximum briglitness using the UV Sclimidt spectrograph. This spectrum 
is classified as "Ia-91bg" by SNID. Clipped bad residuals along with host-galaxy emission lines and telluric absorptions (if present) are 
plotted in light grey. 
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Figure 15. Spectrum of SN 2007ba taken 8 d (rest frame) past maximum brightness using DEIMOS. This spectrum is classified as 
"la-norm" by SNID. Chpped bad residuals along with host-galaxy emission lines and telluric absorptions (if present) are plotted in light 
grey. 




^ I I I I I I I I I 

4000 5000 6000 7000 ^8000 9000 
Rest Wavelength (A) 

Figure 16. Spectrum of SN 2004Y taken on 2004 Feb. 15 using ESI. This spectrum is classified as "la-norm" by SNID. Clipped bad 
residuals along with host-galaxy emission lines and telluric absorptions (if present) are plotted in light grey. 
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Figure 17. Spectral sequence of SN 20071e, classified as "la-norm" by SNID. Each spectrum is labelled with its rest-frame age relative 
to maximum brightness. Clipped bad residuals along with host-galaxy emission lines and telluric absorptions (if present) are plotted 
in light grey. The galaxy-subtracted spectra are overplotted in dark grey. The horizontal dashed lines correspond to zero flux for each 
spectrum. 
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Figure 18. A histogram of the number of spectra versus the number of objects in our SN la sample. Our average is ~2.2 spectra per 
object. 
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Figure 19. A histogram of the redshifts of all of the SNe la in our sample (top) and a zoom-in on those objects with z 0.1 {bottom). 
The shaded regions represent objects for which we have phase information (i.e., a date of maximum brightness). The dotted vertical line 
(z = 0.015) is approximately the redshift above which peculiar velocities can be ignored in cosmological calculations (see, e.g., |Astier| 
|et al.|p006| . The dashed vertical line (z = 0.05) represents our approximate 90 per cent cutoff; i.e., ~90 per cent of our SNe la have 
redshifts less than 0.05. Our average redshift is about 0.0283 and our median uncertainty is 0.00004. 
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Figure 20. A histogram of the phases of all of the spectra in our sample for which we have phase information {top) and a zoom-in on 
those spectra with t ^ 50 d {bottom). Our median uncertainty is 0.38 d. 
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Figure 21. A histogram of the phase of each SN la (for which we have phase information) at the time of our first spectrum (top) and a 
zoom-in on those spectra with tfirst ^ 30 d {bottom). Our median uncertainty is 0.38 d. 
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Figure 22. A histogram of the MLCS2k2 A value (which is a measurement of the light-curve shape) of each SN la (for which we have 
light-curve shape information). Our average A is about 0.12 and our median A and uncertainty is about —0.03 and 0.035. The different 
shadings correspond to each object's final SNID classification (Section |5. 4. 2[ l 
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Figure 23. The host-galaxy redshift versus the MLCS2k2 A value of each SN la (for which we have light-curve shape information). 
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Figure 24. The phase versus the MLCS2k2 A value of each spectrum for which we have light-curve shape information (top) and a 
zoom-in on those spectra with t ^ 90 d (bottom). The median error bar in both directions is shown in the top-right corner of each plot. 
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6.2 Object (Re-) Classification 



6.2.3 SN 1998cm 



Some of the SNe la presented here were originally classified 
as other types of SNe or remained unclassified prior to this 
work. Using our SNID classification scheme we have reclas- 
sified these objects as bona fide SNe la. Similarly, there are 
a few objects in our dataset that, again after applying our 
spectral classification procedure, were found to be examples 
of some of the peculiar SN la subtypes. All of the objects for 
which (re-)classifications were made are described below. 



10, by [Germany et aT] 
"within a week or two 
of maximum brightness" (Germany et al. 111998 1. Our SNID 



This SN was discovered on 1998 Jun. 
d 19981), and classified as a SN la 



classification reveals that SN 1998cm is 91T-like. Figure [27| 
shows this SN compared to its best-matching SNID tem- 
plate, the 91T-like SN 1997br ~8 d past maximum bright- 



6.2.i SN 2000 J 



This SN was discovered on 2000 Feb. 4, by [Puckett et"aL] 
( 2000 1 . Nearly 6 weeks later it was classified as a SN II based 



on the noisy spectrum presented in this work (Filippenko 



& Coil 20001. However, a SNID fit to the same spectrum 



reveals that it is more likely a normal SN la. Figure [28| shows 
this SN compared to its best-matching SNID template, the 
normal Type la SN 1994D ~54 d past maximum brightness. 



6.2.5 SN 200 les 



6.2.1 SN 19910 



This SN was d iscovered on 1991 Mar. 18, by [Mueller fc Fil-| 
ippenko ( 1991| , and classified as a SN la "about 1-2 months 
past maximum brightness" (Mueller & Filippenko 19911. 



Our SNID classification reveals that SN 19910 is 91bg- 
like. Figure [25] shows this SN compared to its best-matching 
SNID template, the 91bg-like SN 2006em ~21 d past maxi- 
mum brightness. 



This SN was discovered on 2001 Oct. 7, by p] ( [200T| ) , but it 
has remained unclassified until now. From a SNID fit to one 
of the spectra presented in this work, we determine that it 
is likely a normal SN la. Figure [29] shows this SN compared 
to its best-matching SNID template, the normal Type la 
SN 2004fz ~22 d past maximum brightness. 



6.2.6 SN 2002bp 



This SN was d iscovered on 2002 Mar. 8, by |Puckett fc Lan-| 
(20021, but it too has remained unclassified until 



goussis 



now. From a SNID fit to our spectrum presented here, we 
determine that it is a SN 2002cx-like object. Upon further in- 
spection, it seems to be a better match to the quite peculiar 
SN 2008ha ( [Foley et aL]|2009a| ) than to th e more "normal " 
members of the SN 2002cx-like class (e.g., [Jha et al.[[2006 |. 
A comparison between our spectrum of SN 2002bp and our 
spectrum of SN 2008ha from 7.5 d past maximum bright- 



ness (previously published by Foley et al. 2009a I is shown 



in Figure |30) We include SNe 2002bp and 2008ha here in 
our SN la sample even though there is uncertainty regard- 
ing whether SN 2008ha was in fact a SN la (e.g. 
et al.||2009| ). 



Valenti 



6.2.7 SN 2004br 



6.2.2 SN 1993aa 



This SN was discovered on 1993 Sep. 19, by [PoUas et al.[ 
(19931, and classified as a SN la "probably about 1 month 
past maximum brightness" ( [Pollas et al.[[T993| . Our SNID 
classification reveals that SN 1993aa is also 91bg-like. Fig- 
ure [26] shows this SN compared to its best-matching SNID 
template, the 91bg-like SN 2007ba ~ 8 d past maximum 
brightness. 



This SN was discovered on 2004 May 15, by [Graham fc Foley[ 
( 2004 1 , and classified as "an unusual type la supernova, simi- 



lar to the spectrum of SN 2000cx" ( Gerardy et al.i2 004). The 
SNID classification of our earliest spectrum of SN 2004br re- 
veals that it is likely 99aa-like. Figure [3T] shows this spectrum 
compared to its best-matching SNID template, the 99aa-like 
SN 2008ds ~5 d before maximum brightness. We are unable 
to confidently determine the subtype of the two older (both 
> 2 weeks past maximum brightness) spectra of SN 2004br 
in our dataset. The best-matching template of one of these 
other spectra is another 99aa-like SN, and the best-matching 
template of the other is a 91T-like SN. Furthermore, Gane- 
shalingam et al. (in preparation) have determined that the 
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Figure 25. Our spectrum of SN 19910 (black) compared to the 91bg-likc SN 2006em (red) ~21 d past maximum briglitness. 
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Figure 26. Our spectrum of SN 1993aa {black) compared to the 91bg-like SN 2007ba (red) ~8 d past maximum briglitness. 
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Figure 27. Our spectrum of SN 1998cm {black) compared to the 91T-likc SN 1997br (red) ^8 d past maximum brigiitness. 
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Figure 28. Our spectrum of the previously misclassified SN 2000J {black) compared to the normal Type la SN 1994D {red) ~54 d past 
maximum brightness. 
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Figure 29. Our spectrum of the previously unclassified SN 2001es (black) compared to the normal Type la SN 2004fz (red) '~22 d past 
maximum brightness. 
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Figure 30. Our spectrum of the previously unclassified SN 2002bp (black) compared to the peculiar (possible) Type la SN 2008ha (red 
7.5 d past maximum brightness, which was previously published by |Foley et al.| ( |2009a[ |. 
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MLCS2k2 A parameter of SN 2004br is -0.152. This spec- 
tral and photometric information increases our confidence 
in the 99aa-like classification of SN 2004br, though it is still 
uncertain. As seen in Figure [31] there are a few obvious dif- 
ferences between its spectrum and that of SN 2008ds; further 
study of this object is warranted. 



6.2.8 SN 2005dh 



This S N was discovered on 2005 Aug. 10, bylMo ore fc Li| 
( 2005[), and class ified as a SN la "near maximum light" 



( Salvo et al. [2005 1. The SNID classifications of both spectra 
of SN 2005dh in our sample reveal it to be 91bg-like. Fig- 
ure [32] shows our earlier spectrum of this SN compared to 
its best-matching SNID template, the 91bg-like SN 2006cs 
~ 2 d past maximum brightness. This object was specifi- 
cally noted to have an "unusually high" expansion velocity 
of 16,000-16,600 km s~^ (based on the minimum of the Si II 
A6355 absorption fea ture; |Salvo et~al.^ 2 005 Gal- Yam et al, 
|2005[). However, both |Salvi3 et al.| ( |2005y and |Gal-Yam et al. 



( 2005 1 used the host-galaxy redshift presented by Falco et al 



1 19991), z = 0.038, as opposed to the actual host-galaxy red- 
shift of 2 = 0.015 ( [Adelman-McCarthy et~aL]|2558| ) . Using 
the correct redshift, we calculate a relatively normal expan- 
sion velocity of ~9300 km s~^ for our spectrum of SN 2005dh 
from a similar epoch. 



All of the spectra referred to in Table [TT] were obtained 
from either Lick or Keck Observatory with the exception of 
the spectrum of the host of SN 2003ah. On 2008 Dec. 28, 
we obtained a 900 s medium-resolution spectrum of the host 
galaxy of SN 2003ah with the MagE spectrograph ( Marshall] 



et al. 2008 1 on the Magellan Clay 6.5 m telescope. Data 



reduction was similar to the process described in Section [3] 
with the exception of sky subtraction. For this spectrum, 
the sky was subtracted from the images using the method 
described by Kelson (20031. Further details of MagE data 



reduction are described in Foley et al. ( 2009a I 



Two other objects in the BSNIP sample which lack pub- 
lished spectroscopic host-galaxy redshifts also deserve spe- 
cial mention. SN 2001ei was discovered in a faint host (LO- 



TOSS J235102 .95-H271050.6) with no known redshift (]Mar- 
]tin fc Lill200T] ); however it is likely within the Abell 2666 
galaxy cluster {cz « 8040 km s~^), possibly associated with 
NGC 7768 [cz « 8190 km s"^). These redshifts are slightly 
lower than the SNID-determined redshifts for our two spec- 
tra of SN 2001ei (as fisted in Table [7|. SN 2004fy was dis- 
covered in MCG +15-1-10, which is probably interacting 
with NGC 3172 {cz ^ 6100 km s~^). This matches weU 
with the SNID-determined redshift for one of our spectra of 
SN 2004fy, though it is somewhat larger than that of the 
other spectrum. Due to the uncertainty of both of these ob- 
jects' host redshifts, they are not listed in Table [l] 



6.2.9 SN 2008Z 



This S N was discovered on 2008 Feb. 7, by [Puckett etliL] 
1 20081, and classified as a Type la (]Blondin fc Calkins|2008l). 



The SNID classification of our earliest spectrum of SN 2008Z 
reveals that it is 99aa-like. Figure [33] shows this spectrum 
compared to its best-matching SNID template, the 99aa-like 
SN 2008ds at maximum brightness. Ganeshalingam et al. (in 
preparation) have determined that the MLCS2k2 A param- 
eter of SN 2008Z is -0.152. 



6.2.10 SN 2008m 



This S N was discovered on 2008 F eb. 13, by [Boles fc Li| 
(20081, and classified as a Type la (Silverman et al. 20081 



using the earliest spectrum of this object presented here. 
Our SNID classification of this same spectrum reveals that 
it is actually 91bg-like. Figure [34] shows this spectrum com- 
pared to its best-matching SNID template, the 91bg-like 
SN 2007ba ~5 d past maximum brightness. 



6.3 New Redshifts for Individual Objects 

Some of the objects in our dataset do not have published 
spectroscopic host-galaxy redshifts. Therefore, we have ob- 
tained host-galaxy spectra of several SNe presented in this 
work in order to determine their redshift. Furthermore, we 
have calculated the host-galaxy redshift of one of these ob- 
jects with no published redshift based on narrow features 
present in our SN spectra. The SNe for which this was done, 
their host galaxies, the redshifts themselves, and basic in- 
formation about the spectrum from which the redshift was 
determined can be found in Table [TTl These redshifts are 
also listed in Table [T] 



7 CONCLUSION 

In this first BSNIP paper we presented a large, homoge- 
neous set of low-redshift [z ^ 0.2) optical spectra of SNe la. 
584 spectra of 199 SNe have well-calibrated light curves 
with measured distance moduli, and many of the spectra 
have had host-galaxy corrections applied. We also presented 
our observing and reduction procedures used during the two 
decades over which we collected these data, as well as our 
"colour matching" method for removing residual galaxy con- 
tamination. Our relative spectrophotometry was shown to 
be extremely accurate for the vast majority of our dataset. 
How the data are currently stored and will eventually be 
made accessible to the astronomical community was also 
discussed. 

In addition, we described the construction of our own 
set of SNID spectral templates as well as our classification 
scheme which utilises these new templates. Using our clas- 
sification procedure we were able to classify for the first 
time (as well as reclassify) a handful of objects as bona 
fide SNe la. Furthermore, we presented classifications of ob- 
jects as members of some of the peculiar SN la subtypes 
that were heretofore assumed to be "normal." In total our 
dataset includes spectra of nearly 90 spectroscopically pecu- 
liar SNe la. We also determined spectroscopic host-galaxy 
redshifts of some objects where these values were previously 
unknown. 

The sheer size of the BSNIP sample and the consis- 
tency of our observation and reduction methods makes this 
sample unique among all other published SN la datasets. 
In future BSNIP papers we will use these data to exam- 
ine the relationships between spectroscopic characteristics 
and other observables (such as photometric and host-galaxy 
properties) . 
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Figure 31. Our spectrum of SN 2004br (black) compared to the 99aa-likc SN 2008ds (red) ^5 d before maximum briglitness. 
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Figure 32. Our spectrum of SN 2005dh (black) compared to the 91bg-like SN 2006cs (red) ~2 d past maximum brightness. 
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Figure 33. Our spectrum of SN 2008Z (black) compared to the 99aa-like SN 2008ds (red) at maximum brightness. 
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Figure 34. Our spectrum of SN 2008ai (black) compared to the 91bg-like SN 2007ba (red) ~5 d past maximum brightness. 
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Table 11. Previously Unpublished Spectroscopic Host-Galaxy Redshifts 



SN Name 


Host 


'^^hclio 


UT Date 


SN/GaF 


Abs/Emis'^ 




Galaxy 


(km s~^)^ 


of Spectrum'^ 






SN 2003ah 


LOTOSS J044309.01+004553.4 


10153 (3) 


2008-12-28 


Gal 


Emis 


SN 2006mp 


MCG +08-31-29 


8090 (300) 


2006-11-03 


Gal 


Emis 


SN 2008s3° 


2MASX J23004648+0734590 


12300 (300) 


2008-09-08 


Gal 


Abs 


SN 2008s5f 




9290 (300) 


2008-09-22 


SN 


Emis 



''The redshift uncertainty is in parentheses. 

''UT date of the spectrum from which we determined the redshift. 

"Gal" = Spectrum from which we determined the redshift was of the host galaxy itself; 
"SN" = Spectrum from which we determined the redshift was of the SN but contained narrow 
host-galaxy spectral features. 

"Emis" = Emission features were used to determine the redshift; "Abs" = Absorption 
features were used to determine the redshift. 
''Also known as SNF20080825-006. 
^Also known as SNF20080909-030. 



Our sample is also a preview of coming attractions; 
new large transient searches such as the Palomar Transient 
Factory ( PTF; |Rau et al.||2009| [Law et al.|[2009| ) and Pan- 



STARRS ( [Kaiser et al.||2002[ ) will compile datasets similar 
in size to ours in just a few years. 
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